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ABSTRACT 

The reduction of the discrete tones generated by jet engines is 
essential for jet aircraft to meet present and proposed noise standards. 
Present reduction methods based on acoustic duct treatment are plagued 
by high cost, as well as severe weight, efficiency, and maintenance 
penalties. The discrete tones generated by the blades and vanes 
propagate in the inlet and exhaust duct in the form of spiraling 
acoustic waves, or spinning modes. The reduction of these spinning 
modes by the cancellation effect of the combination of two acoustic 
fields is the thrust of this research. 

The Spinning Mode Synthesizer has provided the means for effective 
study of this noise reduction scheme. Two sets of electrical-acoustical 
transducers located in an equally-spaced circular array simultaneously 
generate a specified spinning mode and the cancelling mode. Analysis 
of the wave equation for the synthesizer has established the optimum 
cancelling array acoustic parameters for maximum sound pressure level 
reduction. The parameter dependence of the frequency ranges of pro- 
pagation of single, specified circumferential modes generated by a 
single array, and of effective cancellation of the modes generated by 
two arrays, has been determined. Substantial sound pressure level 
reduction has been obtained for modes within these limits. 

Analysis of the system has shown that adaptive control is nec- 
essary for the closed-loop maximization of cancellation, A two- 
dimensional, perturbation type, gradient-adaptive controller has been 
developed. Analog simulation has yielded favorable convergence and 
load change response. This research has established a firm basis for 
the development of the cancellation system for fan-duct sources. 
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I. Introduction 

1.1 Background to Problem 

The aviation industry is currently facing technologically demand- 
ing and financially burdensome noise restrictions. The Federal Aviation 
Administration has been empowered by Public Law 90 - 411, "Control and 
Abatement of Aircraft Noise and Sonic Boom n , to establish rules and 
regulations for the control of aircraft noise. Certification noise 
limits for new aircraft have since been established and are shown in 
Figure 1. These limits are for a three point certification test. 
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Certified Takeoff Gross Wt. (LBxlOOO) 
FIGURE 1 (Reference 1) 


FAA Maximum Noise Limits for Certification 


As shown in Figures 2 and 3, current aircraft, with few exceptions, do 
not meet the certification requirements. Future reduction of certifica- 
tion noise limits can be expected, 



FIGURE 2 
(Reference 1) 

Takeoff Noise Levels - Present Aircraft 
(3.5 N. miles @ 1000 ft.) 





EFFECTIVE PERCEIVED NOISE LEVEL (EPNdB) 



FIGURE 3 


(Reference 1} 

Approach Noise Levels - Present Aircraft 
(1.0 N. miles on 3° Glide Slope) 



EFFECTIVE PERCEIVED NOISE LEVEL (EPNdB) 



0 


CERTIFICATED TAKEOFF GROSS WT. (LB x IOOO) 
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The early pure jet engines such as the Pratt § Whitney JT3C (for 
the Boeing 707.' s) or the General Electric CJ 805-3 produced broad band 
noise at high sound pressure levels. These unacceptably high broad 
band noises were reduced a few dB's by recourse to multiple tube jet 
nozzles. After much research, it was established that such pure jets 
were inherently noisy with little chance to substantially reduce noise 
without unacceptable thrust losses. 

Since this jet noise is proportional to anywhere from the 6th to 
the 8th power of the jet velocity the next generation of engines uti- 
lised medium bypass ratio fans. This resulted in much higher airflows 
and reduced jet velocities thereby reducing the jet noise problem. The 
early fan engines such as Pratt § Whitney JT3D and JT4D (for Boeing's 
707 and 720, and Douglas' DC8's) produced high levels of noise resulting 
from the high speed fans. This fan noise is characterized by pure tones 
at the blade passage frequency and harmonics giving the typical "whine" 
of these aircraft. 

Extensive research has been conducted on noise reduction schemes 
for these engines. Substantial noise reductions have been obtained by 
proper selection of the primary cycle variables. The design of the by- 
pass ratio and the fan pressure ratio consistent with mission and noise 
considerations has later resulted in quieter engines due to reduced flow 
noise. 

Jet noise and pure tones from the fan, compressor, and turbine 
are of concern. The reduction of the dominant pure tone noises for these 
later engines is especially important for the EPNdB rating where pure 
tones are penalized. The installation of acoustical treatment in the 



7 


engine ducts and acoustically treated splitters has been shown to provide 
effective pure tone reduction, but with high cost, as well as weight, 
efficiency, and maintenance penalties. Current and proposed acoustical 
treatment is illustrated in Figure 4. 

To reduce these adverse effects the possibility of acoustically can- 
celling the pure tones appeared very attractive, provided this feat could 
be accomplished economically. The cancellation device would only operate 
on take-off and landing and thus not effect mission performance except for 


SLS THRUST = 97,000 LB 
BYPASS RATIO = 6 



III! TREATMENT LEVEL 

FIGURE 4 (Ref. 3) 

Current and Proposed Acoustical Treatment 


the additional weight which would conceivably be less that the added 
weight of the acoustically treated ducts. 
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1.2 Spinning Mode Synthesizer 

Pure tone noise reduction schemes have been tested on fan-duct 
systems. Broad band noise and the wide spectrum of modes of sound 
generated in a fan-duct system have made meaningful parameter varia- 
tion studies difficult. The synthesis of these pure tones, or spinning 
modes of sound, would allow the isolation of specific modes. Effective 
analysis of noise reduction schemes could then be conducted. 

In order to accomplish this isolation of individual duct modes, 
the Spinning Mode Synthesizer (SMS) was developed at the Noise Control 
Laboratory of The Pennsylvania State University to generate spinning 
modes of sound of controlled complexity without the generation of 
broad band noise (8, 14). The SMS consists of a circular array of 
electrical-acoustical transducers (commercial loudspeakers) mounted 
on a flat plate, as shown in Figure 5. 


r 



FIGURE S 


Spinning Mode Synthesizer 
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The transducer locations can be specified by: 
z q = axial position of array 

= 0 

r . = radial distance of j** 1 transducer at z = 0 
03 

= r , j * 1, 2, . . . n, where n is the number of 
transducers 

9 . = angular position of transducer at z = 0 
2tt 

= — — (j - 1) radians, j = 1, 2, . . . , n 

The sound, field within the ducts of a jet engine can be described 
by the wave equation for ducts of arbitrary shape. A logical basis for 
the study of the modes of sound associated with such a system would be 
provided by a semi-infinite, constant-radius, rigid-walled duct. This 
acoustic environment has been approximated by placing the Spinning Mode 
Synthesizer at the source end of a terra cotta duct having an anechoic 
termination. This acoustic system, together with the analytical back- 
ground provided by Tyler and Sofrin (4), Lowson (5), and others (6,7), 
provides a basis for meaningful studies of spinning modes of sound. 

A general analytical model has been established for a ducted acous- 
tic system by the solution of the wave equation for a semi-infinite duct 
(See Appendix 1) . For the system shown in Figure 6, 
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r 



FIGURE 6 

Cylindrical Coordinate System 


we have 


V 2 $ (r,6, z,t) 


3 2 » (r>6,z,t) 
9t 2 


where $ is the acoustic velocity potential (Reference 15) and the V 2 
operator for the cylindrical coordinate system is given by: 


y2 _ r 9 / *1 


1 ^ 


r L r2 802 


3z‘ 


An assumed solution of the form 


<J> = R (r) 0 (6) Z (z) T(t) 

with an assumed sinusoidal input, the homogeneous boundary conditions 


1) 0 (0) = 0 (0 + 2t r) 
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2 ) 


30 (6) 
90 


90 (6) + 2tt) 
90 


3 ) 


3R(r) 

9r 


1 r=b 


= o, 


the inhomogeneous boundary condition 


v -loot 

-z = f(r,0,O) e , 

z=0 

and the boundedness condition 
1) R(0) = finite 

yields an acoustic pressure field 

P(r,0,z,t) = i u) P $ 


P 


00 oo 

Z l 
m=0 y=0 


E 

my 



imB 

e 


+ 


H 

TTl]i 


-ime, . 


J (k r) e^^my 2, 
m v my J M 


The amplitude coefficients E , G , and H are given in Appendix 1, 

m y my' my ® 

See nomenclature or Appendix 1 for parameter definitions. 

The circumferential mode number m and the radial mode number y 
are two basic descriptive parameters. The concepts associated with 
these parameters can be best explained by an examination of the appro- 
priate components of the pressure solution. 

The circumferential dependence of pressure can be shown by 


where 


n im0 -im0 

P = g + © + g e 

m is a non-negative integer. 

Re [P] ■+ | g | cos m 9 


This reduces to 


Therefore, the circumferential distribution consists of m cycles of 
a cosine wave, as illustrated by Figure 7. 
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The radial distribution is governed by the boundary condition: 



where R = radial component of solution to wave equation 
r = radial parameter 
b = duct radius 

From the solution of the wave equation for the system, we have 

t-[J (k r) ] , = 0 

Sr 1 m my r=b 


or 


J m < k my b) = 0 


where J = Bessel function of the first kind of order m 
m 


th 

k = eigenvalue of order y for the m mode having a hub-tip ratio 
0 = 0 


The values of k^b, for m and y, that satisfy the equation are given 
in Table 1. 


TABLE 1 


t 

Roots of J (k b) = 0 
m my 


k b 
m]s 

o 

ii 

P- 

i 

y = i 

y = 2 

y = 3 

y = 4 

m - 0 

0 

3.830 

7.012 

10.168 

13.317 

m = 1 

1.840 

5.328 

8.522 

11.700 

14.856 

m = 2 

3.053 

6.703 

9.964 

13.164 

16.339 

m = 3 

4.199 

8.011 

11.340 

14.578 

17.780 

m = 4 

5.315 

9.278 

12.675 

15.955 

19.186 
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The nature of propagation is determined from the term e ' M in 
the solution to the wave equation (see Appendix 1). The wave number, 
in the axial direction, is given by 



For v to be real, 
'my 



or 


f > 


K c 
my 


— 2t r 


The frequency at which Y m y = 0 is defined as the "cutoff" frequency, 
*my ' 


where 


b = 


f 

my 

.5 ft. 


(k b) c 
^ my _ 

2TTb 


and k 


my 


is given in Table 1. 


The "cutoff" fre- 


quencies for m and y are given in Table 2. 
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TABLE 2 

Cutoff Frequencies for One Foot Diameter Duct 


f 

mu 

y = 0 

y = 1 


y = 2 

hO 

II 

P- 

y = 4 

m = 0 

0 

1378 


2524 

3660 

4794 

m = 1 

662 

1918 


3067 

4212 

5348 

m = 2 

1099 

2413 


3587 

4739 

5882 

m = 3 

1511 

2883 


4082 

5248 

6400 

m = 4 

1913 

3340 


4563 

5744 

6906 

Duct resonance occurs at y 

my 

= 0 

or f = f 

my 

. For y 

my 

real, or 

£ ± f mr 

, the exponential term 

i is 

complex and 

the mode propagates. 

For *mu 

comp lex, or 

f < V 

the mode decays 

exponentially with 


axial distance since the exponential term is real and negative. 


The decay rate decreases as frequency increases and approaches 

zero as the frequency approaches the cutoff frequency. 

The radial distribution of the pressure field inside a duct 

with no center section is given by J (k r) where m and y are 

° m my 

specified for a particular mode. The distributions for low order 
m and y are shown in Figure 8. The radial mode number y thus 

p 

this corresponds to the number of times ^ crosses the r - axis 

max 


in Figure 8. 

High intensity spinning modes with circumferential distribution 
m and radial distribution y are generated by assigning all transducers 
equal source strength Q, by operating above the cutoff frequency for 



16 


the specified mode, and by assigning the transducer the acoustic 
phase 

(j - 1 ) radians, j * 1, 2, .... n 
where n is the number of transducers. 



FIGURE 8 

Radial Distributions of Sound Pressure 


A set of modes {M} is defined to contain the modes which can be 
excited by a given source amplitude and phase configuration. Although 
m must be positive for the wave equation analysis, both positive and 
negative m will be allowed in this set. The negative sign will desig- 
nate spin in the opposite direction of the synthesized shaft rotation. 
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Amplitude measurements were made in terms of sound pressure level, 
which is given by 

SPL - 20 log 10 £— 

where p = rms sound pressure 

2 

p ref = .0002 dynes/cm . 

The radial distribution of sound pressure level for low order m and 
y are shown in Figure 9. 


SPL 


m| = 0 



m| =2 



0 



rwYi 

OD ID 



FIGURE 9 


Radial Distributions of Sound Pressure Level 
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1.3 The Cancellation Concept 

The high intensity discrete tones generated by rotor-stator 
interaction propagate in the ducts of jet engines in the form of com- 
plex spinning waves. But, only a few modes associated with any given 
interraction propagate unattenuated. The reduction of these modes 
could therefore provide a substantial reduction of the pure tone noise. 

A method for the reduction of these spinning modes is the subject 
of this research project at The Pennsylvania State University. It has 
been proposed to reduce the sound pressure level of the unattenuated 
spinning modes of sound generated by a jet engine by the cancellation 
effect of the combination of two acoustic fields having equal ampli- 
tudes but opposite phase. Cancellation would be maximized by the 
optimization of an amplitude parameter and a phase parameter of 
another spinning mode source. 

The Spinning Mode Synthesizer provides the means for effective 
study of this noise reduction scheme. In addition, a modified SMS can 
be used as the source for the cancelling acoustic field in actual appli- 
cation, The cancellation of simulated pure tone fan noise forms the 
basis for this investigation. 

The noise source and cancelling source can be simulated simultane- 
ously by the SMS by assigning half its transducers to each source. To 
obtain the necessary symmetry for proper mode excitation, every other 
transducer is assigned to the source array with the interstitial trans- 
ducers being assigned to the cancelling array. The letter "n” shall 
designate the number of transducers in each array (4 for the case of 
an 8- transducer SMS) and the subscript j shall represent the j^ 1 
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transducer of the array in question. The acoustic phase formula 
presented in Section 1.2 applies to the source array. The optimal 
parameter solution for the cancelling array is to be determined. 

By simple waveform analysis, the optimal solution for the 
minimization of the spinning mode for a field generated by a 
circular transducer array source given by 


r . 

oi 


til 

= radial distance of j transducer 


Cl. 3-1} 


r Q ,j=l,2, ...,n 


0 . 
oi 


til 

= angular position of j transducer 


(1.3.2) 


~ (j -1), j = 1, 2, . 




having 


Q o j = source strength of j transducer 


(1.3-3) 


- Q q ,-1,2, . . . , n 


4> o j = acoustic phase of transducer 


(1.3-4) 


2nm , . . . . . . _ 

= — (l - 1) , j = 1, 2, . . . , n 

by a circular cancelling transducer array positioned at 


r.=r ,j=l,2, ...,n 
oj o ’ J * * * 

9 . = — (2j - 1) + - , j = 1 , 2, . . . , n 

oj n J 1 n * J ’ ’ * 


Cl. 3-5) 
(1.3-6) 


is expected to be 


q o j = source strength of j 1 "^ 1 transducer 


(1.3-7) 



> j = 1 > 2 , . . .,n 
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a Q j = acoustic phase of j ^ transducer (1.3 8) 

= (2j - 1) + ir , j = 1, 2, . . n. 

It is assumed that both arrays are driven at the same frequency. Suc- 
cessful test cases utilizing the principles of the above solution were 
obtained previous to this author's investigation. 

1.4 Automatic Control of Cancellation 

The radial and circumferential distribution of spinning waves 
generated by a jet engine depend upon engine parameters and operating 
conditions. The numbers of blades and vanes in a rotor-strator set 
and the engine rpm are related directly to the blade passage frequency, 
f = n’ BN 

where 

f = blade passage frequency 

i 

n = harmonic number 

B = number of blades 

N = angular speed of shaft in cps 

which determines the modes generated by rotor-stator interaction. The 
rotor speed alone determines the modes generated by the rotor only. 
Throttle setting and air speed are major factors in variation of blade 
loading, and thus the pressure fields which generate spinning modes. 

The magnitude of the pressure fields along with duct geometry determine 
the amplitudes of the spinning modes generated. 
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Engine RPM, throttle setting, and air speed are time variables. 

In control systems context, these parameters arc slowly time varying. 
Thus there will be time variation of the dominant modal distributions 
and wave parameters. The reduction of the sound pressure level of 
the dominant modes by a cancellation technique will therefore require 
automatic control . 

The development of automatic control for the Spinning Mode Synthe- 
sizer used in this research will provide the background necessary for 
the application of these approaches to full scale. The automatic con- 
trol design objective will provide an organized approach to the investi- 
gation of the system. The implementation of automatic control would 
enable efficient parameter variation, and thus rigorous evaluation of 
the synthesizer and of the cancellation concept. The design and 
development of an automatic controller is thus the intended scope of 
this research. 

II. Experimental Verification of In-Duct Sound Pressure Level Reduction 
by Cancellation 
2.1 Description of Apparatus 

The Spinning Mode Synthesizer was operational at the start of 
this research. The transducer-duct-measurement system is shown 
schematically in Figure 10. 
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FIGURE 10 

Transducer-Duct-Measurement System 


The duct consists of four terra cotta pipe sections (Figure 11) , 
which provided an inexpensive, acoustically rigid duct. The anechoic 
termination (Figure 12) reduces acoustic boundary effects and provides 
an approximation to an infinite duct with all waves travelling in the 
positive axial direction. The duct diameter is 12 inches. 

The Spinning Mode Synthesizer as developed by John M. Seiner (8) 
consists of the end plate, transducers, and control electronics. 

Figure 13 shows the end plate and transducers. The transducer ampli- 
tude and phase controller is shown in Figure 14. 




Figure 11 


Duct, SMS, and Probe Mechanism 
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Figure 12 

Anechoic Termination 



Figure 13 

Spinning Mode Synthesizer 




Figure 14 

Amplitude and Phase Controller and Instrument Console 
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The basic duct instrumentation is shown in Figure 15. A well 
equipped console was available (Figure 14) so other instrumentation 
schemes were utilized when needed. 

In-duct measurements were made by several systems. The probe 
mechanism (Figure 16) is a microphone probe capable of continuously 
traversing the duct radially from the center to the wall of the duct 
at selected circumferential and axial positions. The abcissa for 
radial X-Y recorder plot distributions was obtained by generating 
an electrical signal proportional to radial distance by a potentio- 
meter circuit (Figure IS) . Two probe mechanisms and two recorders 
were used to accelerate testing. 

Measurements near the transducers were obtained by a plate- 
supported microphone mechanism (Figure 17) . The microphone was posi- 
tioned manually to a point centered on the transducer and flush with 
the protective screen, as illustrated in Figure 17. This results in 
approximately a quarter inch distance from the transducer membrane to 
the microphone diaphram. For the frequency range under consideration, 
500 Hz. to 6000 Hz., the wavelengths exceed the 1/2 inch transducer 
radius, so such a close measurement position is allowed. To position 
the microphone circumferentially the transducer under concern was 
operated at low voltage. The microphone was then positioned by maxi- 
mizing the amplitude of the measurement signal, which was observed on 
an oscilloscope. It is important to note at this time that all sound 
pressure level measurements were made with respect to 1 volt. 
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Figure 15 


Basic Duct Instrumentation 
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Figure 16 
Probe Mechanism 
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Figure 17 

Plate Microphone Mechanism 
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2.2 Experimental Results 

Preliminary experimentation was of an exploratory nature. Suc- 
cessful cancellation test cases had been obtained previous to this 
author's investigation. The experimental determination of the set of 

modes {M }, m = 0, +1, +2, . . y« 0, 1, 2, . . . , which could be 
my 

generated and cancelled by a 4-transducer array was considered basic 
to this initial investigation. It was also desired to establish the 
best open loop cancellation procedure. 

The assumption of an ideal electrical-acoustical system provides 
a logical basis for initial investigation. This encompasses the follow- 
ing specific assumptions: 

1) perfectly round, rigid-walled duct 

2) an ideal anechoic termination so that .traveling waves 
propagate only in the positive axial direction 

3) symmetric transducer location on the end plate 

4) identical transducers 

These assumptions enable application of equations 1.3-3 through 1.3-8, 
governing the acoustic strength and phase of transducers for generation 
and cancellation of spinning modes, to the electrical inputs to the 
transducers. The appropriate form is: 

V . = input voltage to source transducer (1.3-3a) 


= V , j = 1, 2, . . . , n 

= electrical phase to j ^ source transducer 


2 ™ 


Cj - 1), j = 1, 2, 


(1 . 3-4a) 
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v , 
OJ 


input voltage to j 


.th 


cancelling transducer 


(1 . 3-7a) 


= V, j = 1, 2, . . . . , n 
0 o ^. = electrical phase to 


cancelling transducer (1.3-8a) 


= ^ (2j -1) + TT, j = 1, 2, . . n 

To verify the predicted optimal solution for the cancelling array, 
the cancelling array can be controlled as a two-dimensional open loop 
system of the form: 

= v, the first variable, j = 1, 2, . . ., n (1.3~7b) 
e oj = (2j -1) + Ct, j = 1, 2, . . n (1 . 3-8b) 


where a is the second variable 

The parameters v and a are then to be controlled to optimize system 
performance. Further adjustment can be obtained by control of the 
cancelling array as an eight dimensional system, with a voltage and 
phase for each transducer. 

Previous to the investigation of spinning modes, plane wave test 
cases were run. Significant sound pressure level (S PL) reduction was 
obtained utilizing the predicted solution: 
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TABLE 3 

Optimal Solution for Plane Wave 


n 

Source 

V 

Array 

0 

Cancelling 

V 

Array 

e 

i 

1 . 0 volt 

0 rad 

1.0 volts 

tt rad 

2 

1.0 

0 

1.0 

TT 

3 

1.0 

0 

1.0 

TT 

4 

1.0 

0 

1.0 

TT 


A test case with frequency of 1400 Hz is shown in Figure 18. 

The fundamental circumferential spinning mode, Mj p = 0, 1, 2, 

. . was investigated next. The source transducers were assigned 

7T 3 

the relative electrical phases of 0, it, and -tt radians for the 

respective transducers according to equation 1.3-4a. The cancelling 

. . tt 3tt 5tt , 7tt 

transducers were set respectively to j + a, + a, + a, and j- + a 

radians, where 0 <_ a < 2tt. From the idealized far field analysis, the 
cancellation concept, a ^ is expected to be tt radians (equation 1.3-8a). 

All transducer voltages were made equal according to equations l,3-3a 
and 1.3-7a. 

The Mjj spinning mode (m = 1, p = 1, 662 Hz f < 1918 Hz) was 
thoroughly investigated. The desired SPL reduction was obtained for 
frequencies below 1400 Hz. The manipulation of the cancelling array 
as an 8-dimensional system after 2-dimensional cancellation had been 
optimized resulted in almost complete cancellation (Figures 19 and 20). . 
Above 1400 Hz, the system behavior departed significantly from ideal. 

As illustrated by Figure 21, a = tt radians did not provide maximum 
cancellation. Acceptable SPL reduction could not be obtained with any 
a. Investigation of the mode (m = 1, p = 2, 1918 Hz <_ f < 3067 Hz) 



FIGURE 18 


SPL versus Radial Distance at 6’ from Plate 
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FIGURE 19 


SPL versus Radial Distance at 6' from Plate 
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FIGURE 20 


SPL versus Radial Distance at 6' from Plate 


for m = 1 , u = 0 Mode at 


1240 Hz 
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FIGURE 21 

SPL versus Radial Distance at 6’ from Plate 
for m = 1, v = 0 Mode at 1560 Hz 
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did not clarify the undesirable high frequency cancellation perform- 
ance of the system because of the complexity and variation of the 
results (Figure 22). Investigation of higher order radial distribu- 
tions (larger y) resulted in similar unpredictable response. Further 
testing did not contribute to interpretation of the results. 

Manipulation of the cancelling array as an 8-dimensional system 
to improve system performance resulted in an additional complication 
of a "point cancellation" effect. The sound pressure level could be 
reduced to the ambient level at any single point in the duct. But, 
this often resulted in degradation of the overall system performance. 
Thus, the SPL at a single point could not be used as a system sound 
pressure level index of performance. By a trial and error technique, 
the SPL at single points was manipulated until an overall improvement 
in the radial distribution of the SPL was obtained. The testing pro- 
cess would have been greatly accelerated by a measurement system yield- 
ing a signal proportional to overall cancellation performance. 

The limited frequency range of acceptable sound pressure level 
reductions was not considered adequate cancellation performance. An 
extensive investigation of the system was begun to uncover the source 
of this behavior and to expand the knowledge of the system. 

The mass loading by the air of an individual transducer would be 
effected by the acoustic outputs of the other transducers as well as 
the duct. In the ideal system, such loading should be equal on all 
transducers. An unsymmetric duct or variation in the transducer 
responses could alter the acoustic field. This would result in per- 
formance degradation. 



FIGURE 22 


SPL versus 


Radial Distance at 6' from Plate 


for m = 


1, ji = 1 mode at 2180 Hz 
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The symmetry of both duct loading and variation of transducer 
response were investigated. The plate probe described in Section 2.1 
was used to set all transducers individually to 0 rad phase and -38 dB 
relative at 1500 Hz. A wide range of electrical voltages and phases 
were required to produce the same acoustic output at the given probe 
position (see Table 4) . This could have been caused by significant 
transducer response variation or unsymmetric duct loading. 

The results of this test indicated a need for an in-depth study. 
Such a study could facilitate system performance improvement, the 
determination of an appropriate cancellation procedure, and the def- 
inition of the plant for synthesis of an automatic controller. Test- 
ing would involve the individual isolation and investigation of system 
parameters. Then progressively elements of the system would be com- 
bined to study the effects of each addition until the complete system 
has been re-established. 

The study of the plate-transducer system in the free field provided 
by a large anechoic chamber enabled sufficient isolation of system 
components for the basic investigation. Tests were repeated with the 
ducted system to evaluate the duct. Acoustic measurements were made 
in the anechoic chamber with an available microphone system. All mea- 
surements were made near the transducers as described for the plate 
microphone in Section 2.1, so effects due to radiation from the back 
side of the plate were negligible. 

Variation of transducer response was investigated by comparing 
the electrical inputs necessary to produce a 0 rad phase, -34 dB 
acoustic output. The range of electrical phases was 15° (Table 5) . 

The range of variation of electrical voltage inputs was .1 volts or 
10% of the total voltage range of 1.0 volts (Table 5) . Therefore, 
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variation of transducer response is significant and contributes to 
performance degradation. 


TABLE 4 

Variation of Transducer Response 



Acoustic 

Parameters 

Electrical 

Parameters 

Transducer 

SPL 

o 

<t> 

V 

e 

1 

-38dB 

0 rad. 

1.125 volts 

36 rad. 

2 

-38 

0 

.98 

29 

3 

-38 

0 

1.11 

33 

4 

-38 

0 

1.10 

33 

5 

-38 

0 

1.10 

49 

6 

-38 

0 

1.11 

53 

7 

-38 

0 

1.09 

27 

8 

-38 

0 

1.10 
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Transducers were operated individually at 1500 Hz. 


The symmetry of free field radiation was investigated by the 
generation of low order modes by the source array alone, by the 
cancelling array alone, and by the combined' arrays. Near perfect 
symmetry was observed for a plane wave for the spinning wave, and 
for both arrays set for a cancelled spinning wave (Tables 5, 6, 
and 7). Therefore, as expected, free field radiation and thus trans- 
ducer interaction is symmetric. 



TABLE 5 


Symmetry of Free Field Radiation in 
Transducer Near Fields for Plane Wave 


Transducer 

Electrical 

V 

Parameters 

6 

Plane Wave 

Total Array 
dB <|> 

Source Array 
dB <P 

Cancelling Array 
dB <J) 

i 

.628 volts 

-170° 

-34.15 

3.0° 

-33.7 -3.0° 



2 

.660 

-175 

-34.2 

4.0 


-33.68 

-4.0° 

3 

.690 

-171 

-33.96 

4.5 

' -33.62 -2.5 



4 

.562 

-165 

-34.18 

3.0 


-33.78 

-3.5 

5 

.582 

-170 

-34.22 

4.5 

-33.72 -4.0 



6 

.590 

-161 

-34.1 

3.5 


-33.68 

-2.5 

7 

.605 

.■ -172 

-34.1 

4.0 

-33.68 -4.0 



8 

.613 

-170 

-34.2 

4.0 


-33.80 

-4.0 


Transducers were set individually to -34dB and 0° phase at 1500 Hz. 
Measurements were then taken in the transducer near fields while plane 
waves were generated. 
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TABLE 6 


Symmetry of Free Field Radiation in Transducer Near Fields 
for First Order Spinning Mode 


Transducer 
■ n 

Electrical 

V 

Parameters 

0 

Spinning Modes 

Total Array 
AdB A(f> 

Source Array 
AdB A<f> 

Cancelling Array 
AdB A<£ 

1 

0° 

.650 volts 

-175° 

-1.3 

2.0° 

0.0 

-.5 


2 

45 

.685 

-130 

-1.53 

1.5 



-.1 -.5 

3 

90 

.720 

-80 

-1.4 

2.0 

-.02 

-.5 


4 

135 

.570 

-30 

-1.3 

2.5 



0.0 -.5 

5 

180 

.590 

14 

-1.58 

2.0 

-.11 

-1.0 


6= 

225 

.610 

65 

-1.15 

.5 



.07 -1.5 

7 

270 

.640 

97 

• -1.2 

1.0 

0.0 

-2.0 


8 

315 

.640 

144 

-1.3 

1.0 



-.1 -1.5 


Transducers were set individually to -34 dB and the appropriate acoustic phase for 
an m = 1 spinning mode at 1500 Hz. Measurements were then taken in the transducer 
near fields while spinning modes were generated. 



TABLE 7 


Symmetry of Free Field Radiation in Transducer 
Near Fields for a Cancelled First Order Spinning Mode 


Transducer 


Cancelled Spinning Mode 
AdB A4> 

1 

0 

1.15 

-4.0° 

2 

225 

1,05 

-4.0 

3 

90 

1.1 

-4.0 

4 

315 

1.1 

-4.5 

5 

180 

1.1 

-4.0 

6 

45 

1.1 

-4.0 

7 

270 

1.1 

-5.0 

8 

135 

.95 

-4.5 


Transducers were set individually to -34 dB and the appropriate 
acoustic phase for a cancelled m = 1 spinning mode at 1500 Hz. 
Measurements were then taken in the transducer near fields while 
the cancelled mode was generated. 
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The frequency dependence of the variation of transducer response 
was investigated next. Each transducer was individually set to -36 dB 
at 1000 Hz. The frequency was then swept from 500 Hz. to 6000 Hz. The 
input voltage and phase were governed by a compressor circuit, a special 
purpose voltage regulator, greatly simplifying testing and still enabl- 
ing meaningful conclusions. There was appreciable variation between 
transducers over the frequency range for both SPL (Figure 23) and 
phase (Figure 24) . Thus the variation of transducer response from 
unit to unit is also frequency dependent. 

The duct loading on the transducers can be evaluated by a com- 
parison of the free field frequency response with the ducted frequency 
response. The magnitude of duct interaction is expected to greatly 
exceed transducer variation. So, only a single transducer need be 
studied in the free field. 

Figure 25 shows SPL versus voltage for a variety of frequencies. 

An empirical relation can be established from this figure: 

SPL (v,f) = k (f) + k (f) [20 log ^ ] (2.2-1) 

ref 

where k^(f) and k^Cf) are frequency dependent parameters, k ^ ( f ) is 

the SPL as a function of frequency for the reference voltage. 

Figure 26 shows k^(f) for three references candidates. The functional 

complexity of k^(f) can be an important factor in the synthesis of an 

automatic controller. The parameter is the s ^°P e the lines 

for SPL versus input voltage (Figure 25) . is a constant for any 



FIGURE 23 

Sound Pressure Level Versus Frequency of 
Transducers in a Free Field at Constant Voltage 
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FIGURE 24 

Phase Response versus Frequency of 

Transducers in a Free Field at Constant Voltage 

at (0 = 0 . , r = r . , z = O') 

03 03 









FIGURE 25 


Input Voltage versus SPL for a 
Single Transducer in a Free Field 
at ( r = r Q1 , 9 =9^, z = 0) 




FIGURE 26 

SPL versus Frequency for a Transducer in 
a Free Field for Constant Voltages 
at (r = r 0lJ 0 = 0 O1 , z = 0) 
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given frequency so it can be determined from Figure 26 by calculating 
the rate of change of SPL between any two constant voltage curves. 

The results are given in Figure 27. 

The plate was remounted to the duct and the same tests as those 
conducted in the large anechoic chamber were repeated. Figures 28 
through 35 show SPL and phase response, measured near the transducers 
by the plate microphone mechanism, versus frequency at constant voltage 
from 500 Hz to 1000 Hz for all eight transducers. A wide variation of 
frequency response of opposing transducers (1 and 5, 2 and 6, 3 and 7, 

4 and 8) was nearly identical. Re-examination of the free field results 
(Figures 23 and 24) yielded few indications of this behavior. So, the 
duct was decidedly responsible. 

It is important to note the lack of permanent phase change as the 
frequency is swept through the m = 1, p = 0 resonance at 670 Hz. There 
is a clear 90° phase shift as resonance is reached, but the phase 
returns to its previous value after resonance is passed. Examination 
of the SPL response clarifies this apparent dilemma. The SPL and 
phase at the microphone are the result of the transducer parameters 
and the combination of all circumferential and radial modes which can 
be excited in a semi-infinite duct. Decaying modes, as well as pro- 
pagating modes, contribute to the SPL and phase because the measure- 
ment point is near z = 0, so axial decay is small. The SPL rapidly 
drops off after the m = 1, p = 0 resonance, indicating a rapid 
increase in damping of this mode. The corresponding rapid decay of 



FIGURE 27 

k 9 (f) versus Frequency 






FIGURE 28 

SPL and Phase Response versus Frequency at 
Constant Voltage (.9 volts) 
for Transducer 1 of the Ducted Array 
at (r = 9 = z = 0) 


\ 





FIGURE 29 

SPL and Phase Response versus Frequency 
at Constant Voltage (.9 volts) for 
Transducer 2 of the Ducted Array 
at ( r — ^Q 2 > ® _ ®Q 2 * 2 - 0 ) 






FIGURE 30 

SPL and Phase Response versus Frequency 
at Constant Voltage C-9 volts) for 
Transducer 3 of the Ducted Array 
at (r = r Q3 , 6 = 0^, z = 0) 
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FIGURE 31 

SPL and Phase Response versus Frequency at 
Constant Voltage (.9 volts) for 
. Transducer 4 of the Ducted Array 

at ( r - r 0 4- 6 ‘ V Z ' 0) 







FIGURE 32 


SPL and Phase Response versus Frequency 
at Constant Voltage (.9 volts) for 
Transducer S of the Ducted Array 
at (r = r 05 , 6 - 6^, z = 0) 







FIGURE 33 

SPL and Phase Response versus Frequency at 
Constant Voltage (.9 volts) for 
Transducer 6 of the Ducted Array 

at Cr = r 0 6’ 6 = W Z = 03 
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FIGURE 34 


SPL and Phase Response versus Frequency 
at Constant Voltage (.9 volts) for 
Transducer 7 of the Ducted Array 
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FIGURE 35 

SPL and Phase Response versus Frequency 
at Constant Voltage (.9 volts) for 
Transducer 8 of the Ducted Array 
at (r = r Qg, 0 = 0 ng , z = 0) 


sc 
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the amplitude coefficient of this mode results in negligable contribu- 
tion to the overall SPL and phase at the measurement point. Thus the 
phase shift still exists; it merely cannot be detected. 

The parameters of duct geometry were investigated to uncover the 
source of this undesirable response. The transducers were mounted on 
the plate within machining tolerances, so the problem was elsewhere. 
There are four possible dimensions of misalignment of the plate with 
respect to the duct centerline. Angular displacements, illustrated 
in Figure 36, 



FIGURE 36 

Angular Misalignment of Plate with Duct Centerline 


were found to be within a half a degree. This alone could not cause 
the observed behavior. The center of the plate could be displaced 
from the duct centerline as illustrated in Figure 37. But d^ and d^ 


were found to be negligable. 
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FIGURE 37 

Displacement of Plate from Duct Centerline 


The terra cotta pipes which form the duct do not have high dimen- 
sional tolerances. An out-of round duct section could easily be present 
and provide the observed unsymmetric acoustic loading. It was found 
that the third duct section from the plate was oval with major and 
minor axis differing by a half inch. This was a likely candidate 
for the observed behavior. 

Excitation of plate vibrational modes by the transducers could 
also contribute to the observed behavior. The plate is bolted at 
three points as shown in Figure 38, This bolting pattern would allow 
vibration in the basic plate modes, A test was conducted to determine 
the effect and magnitude of such vibration on the acoustic field. 

Figures 39 through 42 show radial distributions of SPL and phase for 



FIGURE 38 


Plate Bolting Pattern 


an M n spinning wave of 2000 Hz. Note that at only one position was 
a desired radial distribution observed (Figure 41) . Automotive under- 
coat was then applied to the plate to reduce radiation from the plate. 
Radial distributions at the same probe positions were again recorded 
(Figures 39 through 42). A slightly greater symmetry of the distribu- 
tions at both probe positions was observed, an improvement obtained 
at the expense of the well defined radial distribution mentioned above. 
Though a performance improvement was not achieved, the test served to 
show the sensitivity of the system to slight parameter variation. 

It was shown that several possible causes for unsymmetric trans- 
ducer loading exist. The magnitude of the contribution of each was 
not investigated in this analysis. 



FIGURE 39 

SPL versus Radial Distance at 6' from Plate 
for m = 1, u = 1, Mode at 2000 Hz at a Circumferential 
Angle of 0 radians with Respect to Transducer 1 
of the Source Array 
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FIGURE 40 


Phase versus Radial Distance at 6' from Plate 
for m = 1, y = 1, Mode at 2000 Hz at a Circumferential Angle 
of 0 radians with Respect to Transducer 1 of the Source Array 
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FIGURE 41 

SPL versus Radial Distance at 6’ from Plate 
for m = 1, y = 1, Mode at 2000 Hz at a Circumferential Angle 
of - radians with Respect to Transducer 1 of the Source Array 











FIGURE 42 

Phase versus Radial Distance at 6' from Plate 
for m = I, \i - 1, Mode at 2000 Hz at a Circumferential Angle 
of radians with Respect to Transducer 1 of the Source Array 




(degrees) 





The evaluation of duct loading on the transducers was completed 
by the determination of the ducted transducer frequency response. 

Figure 43 shows that the empirical equation developed for the free 
field frequency response is appropriate: 

SPL (v,f) = kj(f) 4 k 2 Cf) [20 log ] (2.2-1) 

ref 

Some of the curves are nonlinear, making k^(f) a multi-valued function 
of voltage at some frequencies. But the magnitude of variations is 
small so an approximate value would be adequate. The parameter (f) , 

the SPL for v is shown for a single transducer with v re £ = *9 

volts at frequencies from 500 Hz. to 4000 Hz in Figures 28, 44, and 45. 
Peaks can be observed at the duct cutoff frequencies. So, the fre- 
quency response for the ducted transducer is of the same form for a 
transducer in a free field, but reflects the frequency dependent 
behavior of the duct. The complexity of k^(f) along with transducer 
variation can clearly contribute to cancellation difficulties. 

A new cancellation procedure was devised to account for unsymmetric 
transducer loading and the complexity of the transducer frequency 
response. Transducers were individually set by the plate probe micro- 
phone to the appropriate acoustic phases and equal SPL output. The 
procedure was tested with an spinning wave of 750 Hz. Results 
are given in Figure 46 through 51. Reduction of over 20 dB was 
obtained at all probe positions for the cancelling condition. The 
field produced by the cancelling array is 5 to 7 degrees from the 
desired 180° phase displacement from the source array field. Com- 
pensation by the adjustment of the acoustic phases of the cancelling 



FIGURE 43 


Input Voltage versus SPL for a 
Single Transducer in the Ducted Array 





FIGURE 44 


SPL versus Frequency for Transducer 1 
of the Ducted Array at a Constant Voltage (.9 volts) 
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FIGURE 45 


SPL versus Frequency for Transducer 1 


of the Ducted Array at a Constant Voltage (.9 volts) 




FIGURE 46 

SPL versus Radial Distance at 2' from Plate for m : 
Mode at 750 Hz with Acoustic Parameters of 
each Transducer set by use of 
Plate Microphone 


1, V = 1, 








FIGURE 47 


Phase versus Radial Distance at 2' from Plate for 
0 Mode at 750 Hz with Acoustic Parameters of each Transducer 
set by use of Plate Microphone 
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FIGURE 48 

Phase versus Radial Distance at 6' from Plate 

for y = 0 Mode at 750 Hz with Acoustic Parameters of each Transducer 

set by use of Plate Microphone 
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FIGURE 49 


Phase versus Radial Distance at 6' from Plate for m = 1, 

0 Mode at 750 Hz with Acoustic Parameters of each Transducer 
set by use of Plate Microphone 
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FIGURE 50 

SPL versus Radial Distance at 8 T from Plate for m = 1, u = 0 
Mode at 750 Hz with Acoustic Parameters of each 
Transducer set by use of Plate Microphone 


/ '03 
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FIGURE 51 

Phase versus Radial Distance at 8' from Plate for m = 1, y = 0 
Mode at 750 Hz with Acoustic Parameters of 
each Transducer set by use of Plate Microphone 
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array transducers did not produce an overall cancellation performance 
improvement (Figures 46 through 51) . 

Sufficient information had been obtained so preliminary experi- 
mental analysis was concluded. Satisfactory cancellation performance 
was observed to 1400 Hz for the fundamental spinning wave. The basic 
system parameters had been investigated and complications noted. Thus 
adequate knowledge of the system had been established for the synthesis 
of an automatic controller for in-duct cancellation. 
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III. Automatic Control Problem 
3.1 Need for Automatic Control 

The spinning modes of sound generated by a jet engine are depend- 
ent in part on relatively slowly time varying parameters. The set of 
circumferential modes which can be generated depend only on the blade- 
vane numbers. But, the predominant radial mode number associated with 
any given circumferential mode depends on the blade passage frequency, 
and duct diameter. The sound pressure levels associated with a given 
mode will vary, dependent largely upon throttle setting, air speed, 
and air density for any given engine configuration. Since the system 
is subject to time varying parameters, an acoustic cancellation system 
will require automatic control. 

Application of the cancellation scheme to a jet engine poses a 
further complication. Because of the high acoustic amplitudes associ- 
ated with fan noise, the linear acoustic approximation no longer holds 
and wave steepening occurs as the waves travel down the duct. An 
adaptive control may offer opportunities to compensate for this factor. 

The automatic control of the Spinning Mode Synthesizer is of 
immediate concern. Objectives of this work are the attainment of 
effective parameter variation, rigorous evaluation of the synthesizer, 
and the introduction of complexities characteristic of jet engines. 

Open loop control of cancellation has been found to be insufficient 
in these areas . 

The determination of the optimum input parameter setting to maxi- 
mize cancellation of a particular spinning mode by open loop control 
is tedious. Open loop control does not allow sufficient test cases to 
ensure effective mode cancellation specifically for high frequencies. 
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The ''point cancellation" effect discussed in Chapter II provides an 
additional barrier to productive experimentation. As suggested, a 
measurement system yielding a signal proportional to overall system 
performance would greatly accelerate the testing process. The syn- 
thesis of such a measurement system is a difficult problem. The 
minimization of "point cancellation" effects on the index of perfor- 
mance will most likely require an empirical design technique. This 
would be impeded by the lengthly experimentation time. Thus the need 
for effective parameter variation is firmly established. 

An important measure of SMS cancellation performance is the range 
of circumferential modes and radial distribution over which effective 
cancellation can be achieved. It is essential to the evaluation of 
the SMS to determine the pertinent parameters and the functional 
relationships which result in undesirable system performance. The 
empirical determination of these relationships would be greatly aided 
by an efficient testing system. 

Open loop control of the cancelling array inhibits introduction 
of complexities characteristic of jet engines. For example, there is 
often significant energy in more than one mode. The cancellation of 
multiple modes increases the dimension of the controller to 
n' = c n 

where n = dimension of the cancellation system for a 
single acoustic mode 


c = number of modes 
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In addition, the spinning modes are contaminated by broad band noise, 
turbulent air flow, a short, thin-walled duct, time varying parameters, 
etc. Manually controlled cancellation under any of these conditions 
would be very difficult. Thus, the development of automatic control 
is essential for the successful development of the cancellation con- 
cept. And, automatic control of the Spinning Mode Synthesizer is a 
basic necessity to the development of an effective in-duct mode 
cancellation device. 

3.2 Statement of Control Problem 

A control system is to be designed and developed for the automa- 
tic control of the Spinning Mode Synthesizer cancellation simulation. 
There are five general variables: the circumferential distribution of 

the spinning wave (m) , the radial distribution of the spinning wave (y) , 
frequency (f) , acoustic pressure amplitude (A), and acoustic phase (cx q ) - 
For a given m, y is a function of frequency so y does not require 
regulation or control. The cancelling array is driven from the same 
frequency oscillator as the source array, maintaining a frequency lock 
on the optimum frequency. The term frequency lock is used to describe 
the cancelling system state of being driven exactly at the frequency 
of the source array, eliminating all beating effects. 

Each transducer has an associated acoustic pressure amplitude A^ 
which will be measured in terms of SPL^ , the sound pressure level near 
the transducer, and an acoustic phase a Q ^ , j = 1,2, 3, 4. The electrical 
inputs, to the transducers are the directly controlled variables. The 
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electrical amplitude and phase will be designated by and 0^. 
Schematically, the plant can be shown as follows: 


Duct Field 





FIGURE 52 
Plant Schematic 


4 

The plant is 10-dimensional: f, m, S (6-, v.). But frequency 

j = l 3 3 

is a variable controlled by a source array. The parameter m would 
be selected by the operator or by a pattern recognition system. So, 
once m is selected, the plant reduces to an 8-dimensional system 
with an 8-dimensional input and single output, or performance indicator. 
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3.3 Determination of System Model 

Cylindrical coordinates for the SMS can be established as 
follows : 


r 



FIGURE 53 

Cylindrical Coordinates of SMS 


The analytical description of the pressure field inside a hard-walled, 
semi-infinite duct due to a general sinusoidal source at z = 0 has 
been developed (Appendix 1) : 


VO w 

P(r,0,z,t) = + m g 0 E my ^my 6 


+ im0 ¥f -1TH0-. T x, v 

+ H e 1 J (k r) 
my J m my 


e 1C W - “ tJ 


where 
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E = am P 1 ^ tude coefficients [Appendix 1) 

G my “ f l" f t r fCr - 6 - 0) J m % r) e ’ lm6 dr d6 ’ " * ° 

% - ! T t r f(r ' 9 ' 0) ( V’ e 1 " 8 dr d0 ‘ " * ° 

G op * H op = C 'o T fCr - 9 ' 0) J o ‘ k op r) dr d6 > " ' ° 

J (k r) = Bessel function of the first kind of order m 
m my 

th 

k = eigenvalue of order y for the m mode having 
my 

a hub -tip ratio a = 0 
f(r,9,0) = amplitude distribution of source 


Y my = / (^) 2 - k 2 
H v c J my 

m = circumferential . distribution 
y = radial distribution 

To complete the solution, the source must be specified. The two 
sets of n transducers will be approximated by two sets of n point 
sources. Both transducer arrays are located a distance r Q from the z 


axis. The relative position of the transducers in the 0 direction 
is given by 


% = ¥ (j " 1} ’ j = 1 > 2 > * • •» n 


. th 


where 6 . is the relative positive of the j transducer. The Dirac 
oj 


delta function is used to mathematically provide the point source des- 
cription in the r and 9 directions. A point in cylindrical 



. 1 14 


coordinates is given by 

g(r,9,z) = i 6 (r - r o ) 6 (8 - 9 q ) 6 (z - z q ) . 

So, a point at r and 0 . in the plane z = 0 is given by 
r 0 0] 

g. (r,9,0) = i 6 (r - r o ) 6 (9 - 9^) 

The cancelling array is offset by an angle A: 


(j-|) source 



FIGURE 54 

Angle Between Cancelling Array and Source Array 


where A = \ [^L (j - D - % (j - 1) - 1) ] = £ . 

The acoustic phase of the source array, using the first transducer as 
a reference, is given by: 

2 ™ 

4> 0 j = (j - 1) , j = 1, 2, . . n 
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th 

where <j> . is the acoustic phase of the j transducer and m is the 
oj r o 

desired circumferential distribution. The acoustic phases of the 
cancelling transducers will be given by a Q ^ , j = 1, 2, . . . , n, 
which are not yet determined. Let Q = source strength of source 
transducer, in cubic meters per second. Let q = source strength 
of cancelling transducer, in cubic meters per second. The source 
can then be defined by: 

f 0,0,0) = Z [£ e 03 6 (r - r Q ) 6(0 - 0 Qj ) 


ia . 


+ r e 03 6 0 - t o ) 6(0 - 9 qJ - J) ] 


Then, 


mu 


^ n ~ i(j) . 

/ 211 / b r { S [— e ° 3 6(r 

° ° j=l r 


r ) 6(e - e .) 
o Oj 


ia 


+ a e “J t(r - rj 6(6 - 9 . - £)1) J(k „r) 


-im9 


oj n' nr my 


dr d9 


- / J (k r) 6(r - r ) dr 
o nr my cr 


I l 2 ; [Q e 1 ^ 6(9 - 9 ) * , 6(6 - 6 - 1) ] e' 1 " 6 de 

j=l 


n 


= J (k r ) Z [Q e 
m my ; o . , 

J = 1 


id) . -im0 ia . -im(0 . + — ) 




oi 'oi n -1 

+ q e J e J J 


Similarly, 


n i<j> . im6 . ia . im (6 . + — ) 

H = J (k r ) .1 [Q e ° 3 e 0:1 + q e 3 e 3 n ] 

mu nr my o j = l ^ 
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Also, 


G + H = / 2lT f h t {.£, [- e' 6(r ~ r ) 6(6 - 0 ,) 
opt opt oo j = l l r o o} 


ia . 


+ a e °i 6(r - r 0 ) 6(6 - 6 oj - J) ] > J 0 ( k oy r) dr d0 


f o J o (k ou r) 6(r ‘ r o ) dT * 


n id) . ict * ^ 

f 2 J [Q e 0:1 6(6 - 0 Qj ) + q e ° 3 6 (6 -6 oj - -) ] d0 


n it)) 


ia . 


■ V k ov, r) jil [Q e 03 ‘ q e 031 


Now, 


K = G e im0 + H e" im0 
mpi mpt mpi 


= J 


m 


n i (cb . - m0 .) 

<VPi*i {lQ * 03 tqe 


i(a . - m(0 . + —}) 
oj oj n , 


im 0 

e + 


i(4> rt , + mQ n .) ' i(a n . + m(S Q . + -)) 


[Q e + 


q e 


oj 


] e } 


foT m / 0, For m = 0, 




K * J (k r ) tQ e OJ + q 
my o v oy o' j = l 


ia . 
oj 


e " J ]. 


an expression for valid for all m can be given by: 


n 


K 


my 


- „ i(6 . - m0 J 

= JH J (k r ) .S. { [Q ] 3 + q 

2 m mpt o j = l e 


it-, 


i(a . - m(0 . + — )) 

, oj oj rr -j 


if6 . + m0 .) i(a . + m(0 . + — )) . Q 

e in9 *[qe °3 *qe OJ 03 " ] e- 1 " 9 } 
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where 


e = 

IB 


1 for m = 0 

2 for ib 1 . 

Thus the complete solution is given by 


P(r, 9 ,z,t) + K my J m^ k my r ^ e 


i(a)t - 


y z) 
'my 


where K is given above. Since (j> . = m 0 . , K can be written as: 
my oj o oj m y 

m 7 T 

p n ifa . - <& . - 1 ini 9 

K ( k r) .S. {[Q + qe ° J " ]e ° 

my 2 m v m y o J j = l LH J 

o oo 


m it 

2i6 . ifa . -*■ <t> . + —1 -m 8 

■ Y nn v m T m nJ O 

[Q e J + q e 


°j ’°J n 7 je ° ) 


Maximum cancellation occurs when the magnitude of the real compo- 
nent of pressure is minimized with respect to and q, Re[P(r , 0 , z, t) ] 

can be evaluated as follows: 


PCr,e, 2 ,t) * J 0 v So VlK 


n i(<j> . - m 0 .) 

(k r ) .E { [Q e 0J 3 + 

v my o' j=l L 


+ q e 




i(4> . + m0 .) 


bj n" ]e y [Qe -" 0 j -or + 


+ q e 


«■-. * -C» oj ♦ 3» e -« }J (k T) e -,p 

J m v mil J 




iCY^-z - wt) 


m my 


P(r, 6 , z,t) m fo y =0 E my 2 J m^ k my r o'* J m^ k my r - ) 


n 


i(d) . - m 0 . + m 0 + y z - wt) 

' T ai a 'i mil 


• & ^ * 03 05 


my 


i(a . - m 0 . - — + m 0 + y z - wt) 
^ oj oj n my 


q e 


+ Q e 


i f 6 . + m 0 . - m 0 + y z - wt) 
Y oj oj my 
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Then 


i(a . + m0 . + — - m0 + Y™,, 2 - tut) 
oj oj n my j 


+ q e 


m=0 y=0 E my 2 J m^my r o-* J m ^my r ^ 


Re[PCr 1 0 > ,z,t)] = ^ 

I 

* jlj (Q cos (4» o j - m0 Q .. + m0 + Y my 2 - to t) + 


n 


+ q cos (a . - m0 . - — + m9 + Y m ,. z _ + 

M v oj oj n my 


+ Q cos (d> . + m0 . - m0 + y z - tot) + 
H ^ r oj oj my 


+ <* cos (o o1 + m9 oi + ~ - me + y 2 - wt)} 


oj oj 

for m, y such that k £ 


5 ] 


_ r i t 

[m=0 U=0 


U V “5 J m (k myV V k my r) e 


iy z 
'my 


* {Q sin (d> . - m0 . + m0 - tot) + 

j=l oj oj 


mir 


+ q sin (a . - m ® 0 j - — + m0 - tot) + 


+ Q sin (d> . + m0 . - m0 - tot) + 
x VY oj oj 


mil 


+ q sin (a . + m0 . + - m0 - tot)} 

n oj oj n 


to 


for m,. y such that k > — 
’ my c 


!L I n £ J (k r ) J (k r) 
m=0 y=0 m m my o m my 


n 


* [[E my jll {Q cos C % +' W." Wt) C ° S (m6 " m9 oi^ + 


Re[P(r,0,z,t)] 
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. + q cos (a oj + Y mv z - wt) cos (m6 - m0 Q j - ~) > 

foT m, y such that k <_ ~ 
iy z n 

-ME e mU {Q sin C4> 4 - wt) cos ( m0 " m6 m -) + 

L my j = l 03 °J 

+ q sin Cct Q j - wt) cos (m0 - m0 o j - ~) ] 

for m, y such that k > ^ 1 
9 my c 

co , 

considering only propagating waves, ie., k — » we have 

Re[PO,0,z,t) ] = m ? Q ylo E my e m J m^ m y r o^ J m^ m y r ^ 

n 

• {Q cos C4> 0 j + y mv ,z - wt) cos (m0 - m0 Q; .) 

+ q cos (a Q j + y^z - wt) cos (m0 - mO^ - ”)) 

(3.3-A) 

As n -+■ «•, the two arrays approach "source rings", continuous 
arrays in the 0 direction. For this case, Re[P(r,0,z,t)] is zero 

term by term for all m and y independent of r,0,z, and t for 

q = Q and = <J> Q j + (2£.*1) TT, S, - 0, + 1, + 2, . . . . In the 

other extreme, for small n the point source behavior of the arrays 
becomes significant. Re[P] in general cannot be reduced to zero. 

As m -> n, significant multi-mode excitation can make cancellation 
impossible. 

Consider the case of n = 4. As illustrated in Figure 55, there 
is good waveform phasing bias only for m = 0 and m = 1. 


1 
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Acoustic Phases for m = 0, 1, and 2 for n = 4 


For m = 2, spinning waves of equal strength are excited in the clock- 
wise and counter-clockwise directions. Also, the (1,3) and (2,4) 
transducer pairs generate plane waves of equal strength. Thus it 
can be deduced that for a single, well defined waveform, it is 
desirable to have n » m. 

Assuming a single dominant acoustic mode of m = m Q and y = y(f) , 

it is desired to determine the values of q and a . for maximum can- 

cellation. One basic procedure is to take the respective partial 

derivatives and set them equal to zero, thus determining the maxima 

and minima. Re[P] is the real component of a complex number, so both 

positive and negative values must be delt with. The minimization of 
2 

(Re[Pj) provides a physically meaningful and analytically expedient 
means of minimizing the magnitude of P. The sound pressure level. 
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given by 


SPL = 10 log (E ) 2 

Pref 

2 

is minimized by minimizing the rms pressure squared, p . Since 


P 2 = CRe[P]) 2 , 


2 

minimizing (Re[P]) minimizes the sound pressure level. Consider the 
determination of the extremum of {Re [P(m o ,p(f) ,r, 0 , z , t) ] } by evalua- 
tion of the partial derivations with respect to q and a Q _. , j = 1, 

* * 

2, . . . , n. The optimal values q and are then q and 

2 

a Q j evaluated at the minima of {Re[P]}, if and only if these values 


are independent of r,0,z and t. 


2 . 


The partial derivative of {Re[P]} with respect to q is: 


{Re [P]} 2 = 2 Re[P] ^ Re[P] - 0 
£ 

Either Re[P] = 0 or Re[P] = 0. For the acoustic phases, 


fa {Re[P]} 2 = 2 Re[P] |j Re[P] = 0 
oj oj 


If Re[P] = 0, there is only one equation, while there are (n + 1) 

unknowns. Thus if Re[P(m o ,p(f) , r,0,z,t)] is zero, this method of 

* * 

solution cannot be used to determine the optimal values q and . 
If Re[P] i 0, there are (n + 1) equations. 
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sf- - 0 =-> jij (cos (o oj * y - (St) • 

m it 

• cos (m 0 - m l 0 . - ) = 0 

^ o o oj . n 

and 

m 7T 

Jr Re.[P] = 0 =*> cos (a. + Y m z - «t) cos (m Q 6 - m Q 0 Qj - = 

° u oj J o p 

But, the (n + 1) equations give only space and time coordinates 
at which minima can occur in the sound field. Thus no usable infor- 
mation is yet obtainable about q and ct o j . 

To obtain usable information, an index of performance is defined 
as follows: 

ip = . . . 1 i f 2 / 2lT / b {Re[P(m,p)]} 2 r dr d0 dz dwt (3.3-B) 

2. 2_ o o o o 
2J2 ,tt b Z 

This serves as an indicator of the time average pressure amplitude of 
the (m,y) mode over the entire duct volume out to a distance 2. The 
minimization of this index of performance is thus physically meaningful 
since it is the objective of this system to reduce the overall sound 
pressure. Since direct dependence upon spatial coordinates is elimi- 
nated, meaningful results can be expected. The IP is evaluated as 
follows : 

{RetPCm.n )]} 2 - E 2 u e 2 jX/o’ ' 

n 

* [ jSj {Q cos (cf> oj + Y mu - wt) COS (me - mO^) 



123 


2 

+ q cos (a . + y z - wt) cos (m0 - m0 . - —jr) ] 
n oj 'my oj n 


f* { Re[P(m,y)]} 2 r dr = E 2 y e 2 j2 (k^) [•] / q r (k my r) dr 


where [•] represents the bracketed term in the previous equation. 
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i° r J 2 (k r) dr = [ 
J o m my 


1 ^1/ - J l <V. r) * 


2k 


my 


791 2 1 

+ k 2 r (J (k r) ) “] 
my m k my 


(Ref. 10) 


= -Vi<b 2 -, 2 ) jJ ck 2 / w ; c V )) 2 


m 2 J 2 (0)] 
m 


2 k 


my 


2k 


V oi, » 2 - * 2 > <v> 




£ <Re[P(m,p)]) 2 * <* = E 2 u — §- (k^b 2 - » 2 ) J 2 Ck mu b) • 

2 mu 

• Vo> [ -’ 2 

. F 2 t -] 2 
my L J 


where 


?l, - i CE„, -t=— ) 2 (k 2 ,b 2 - n, 2 ) J 2 (k b) J 2 Ck mp r o ) 


my 2 v my k my ' '"my 


Next, 


/f /J {Re [P(m.p)]} 2 rdrdQ - F 2 p /f [-] 2 de 
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n 


F 2 / 2Tr [ .E, (Q cos (<f> . + y z - cot) cos (m6 - m 0 .) 
my o L oj 'my °J 

+ q cos (a . + y ,,z - cot) cos (m0 - m0 . - —7} ] 2 d0 
n oj my c 


oj n 


- C [ k ll jh '0 “ s (+ok * V,- Ut) COS (me ' me »k’ 


n n 


urn. 


+ q cos (a ok + Y my z - cot) cos (m0 - m0 ok - n 


0 > 


{Q cos (4> oj + Y my z - «t) cos (m0 - mB^) 


ta ni + Y imi z " cos (m0 ‘ m6 oj d0 

z - cot) cos (<fc . + y m .z - cot) 


+ q cos v-pj • > my “ -- ™ oj 

F my ^o *k=l j = l ^ COS ^ok + Y my“ "" VT oj 'my 

• j [cos ( 2 m 0 - m( 6 ok + 6 oj )) + cos m( 0 ok - 0 ^) + 

+ ^ cos ^ ok + Y my Z " Wt) COS ( %j + V ' Ut) ‘ 

• i [cos ( 2 m 0 - m ( 0 Qk + 0 oj ) - ” ) + 

+ cos m (0 ok - e oj )3 + Qq cos (4> ok - Y my z - «t) cos (a Qj ♦ Y my z 

• 7 [cos ( 2 m 0 - m( 0 o . + 0 Qk ) - ^ > cos (m 0 ok - m 0 oj - —)] 


+ Q q cos (♦ oj + Y my z - o>t) cos (a ok + Y my z - • 

* J [cos (2m0 - m (0 Qk + 0^ - — ) + cos (m0 ok - + ^)]])d0 

* J, jlj fir cos m Ce ok - 6 oj ) [Q 2 cos 0f ok * y^z -ct) • 


cot) ♦ 
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• cos «oj *W - wt) * q2 cos (a ok ' “ t) ' 

• cos (a ok ♦ Y mu z ut)]* It Q q cos (i»6 ok - »0 oj - — ) 

' c ° s W 0 k * V ‘ “ t) “ S ( “oj * V " ut > * 

_ , - - mn\ 

+ tt Q q cos (m0 , - m0 . + — J • 

H n ok oj n 

• cos (* 0 . ♦ Y m jj Z - ut) COS (o ok * Y mkl z - tot) } 

Integrating next with respect to (lot) , 

SL2Tf 27Tb 9 

/ f f Re [P] rdr dwt = 


£2 tt 


n n 


/ o W kli jh {cos " (e ok - V ' 

• R 2 \ (cos M ok 1 * oj ♦ 2Y ml , z - 2ut) ♦ cos M ok - * oj )] * 


♦ q 2 i [cos (a ok * a oj . 2 Ylnv z - 2wt) ♦ cos (a ok - a oj )]] * 


nvrn 1 


+ . Qq cos (m0 . - m0 . - — ■) — [cos (d> . + a . + 2y z - 2cot) + 
ok oj n' 2 L ^ T ok oj 'my 1 


oj 'my 


+ cos (<j> , - a .)] + Qq cos (m0 t - m0 . + ~) 
T ok oj' J xn ^ ok on n J 


oj 


oj 


tt [cos (a , + 4> • + 2y z - 2ait) + cos (a . - 4> .)]) du>t 
2 1 ok T oj 'my ' v ok r oj J J 


n n 


F** 2 * ,S, { 

my 


2- .2 {cos m (0 . - 0 .) [Q cos (<J> , - d> .) + 

=1 j=l v ok oj' vr ok y oj' 


+ q cos (a ok - a o .)] * Qq [cos C»0 oj - »6 0j - cos C* ok - a oj + 


+ cos (m0 . ~ th0 . + — ) cos Co i - 4> .)]} 
ok oj n J L ok T oj^ J 
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Since the above expression is independent of z, it is the average 
with respect to z; ~ f 1 [• ] d z = [• ] . Thus the index of performance 

Lj 0 

is : 

IP = ull c ° s " < 0 ok - V K 2 C0S W ok - V + 

* , 2 cos (a ok - o oj )] * Q q [COS (n,e ok - n^. - =-) cos - a Q .) * 

♦ cos (>" 0 O k - m0 O j 4 ir> cos ( %k - V ]) 


To determine the optimal q's and a^^'s, the partials are evaluated 
and set equal to zero. 


k " " F mU ^2 kll jh (2t > [C0S " ( ®ok - 6 0j* “ S ( “ok - “oj 51 * 


n n 


. Q [cos Cme ol - m9 o . - =•) cos (* ok - a oj ) * 


* cos [m6 ok - me * =■ cos Ca ok - * oj )]> 


0 IP = — — [2 F 2 — Kr , L. cos m [9, - 9 . ) 
3a . 9a . L my „,2 k=l ok oi 

Ol 01 Zo ^ 


• [Q 2 cos (4» ok - * oi ) + q cos (a ok - a Qi )] + 


+ Qq I cos Cm6 n . - m6 - — ) cos (A, - a ) + 


ok oi n 


ok oi 


+ cos f m0 ok • m6 oi + COS (a ok ' ^0i )]} + 


+ F my JT {q2 + q2 + 2Qq C0S r COS Ca oi " a oi )}] 
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For 


2 1 r 2 

= 2 F — jr I {q cos m (6 . - 0 .) sin (a , - a .) + 

my 2tj 2 n v ok or ok oj J 


+ Qq .[cos 06 . - m0 . - sin (<f> . - a .)]} 
ok oi n ok oi 1 


n = 4, the expanded form of the partials is as follows: 


7 ^-j- IP «> 0 = 4q { 2 + cos ™ cos (a Q9 - a Q1 ) + 


3-JT 


+ cos mu cos (&Q 3 - a^) + cos m — j cos (a^ 4 - a^) + 


miT , . . 

+ cos — j cos (a^ - + cos nvrr cos (a^ - a^) + 

+ cos _ -J cos C a Q4 - a Q3 ) + 2Q {cos [cos a Q1 + 

+ cos (a 02 - m Q ~) + cos (a 03 - my) + cos (a Q4 - m Q ^y) + 


7T 3tt 

+ cos (a Ql - m Q j) + cos (a Q2 - my) + cos (a Q3 - m Q -j)] + 


+ cos m [cos a^ 2 + cos (a n1 - m^Tr) + cos (a n T _ 111 


01 


03 o 


3tt 

+ cos (a Q2 - rn Q — ■) + cos (a Q4 - my)] + 


Sit , , 3tt. 

+ cos m — ^ (cos a^ 3 + cos (a Q1 - m Q — -) + 


+ cos (a 04 - m o y)J + cos m ~ cos a Q4 > 


9a 


01 


2 7T 

IP ==> 0 = q [cos m y sin (a^ - a^) + cos imr 


3tt 

sin (a Q3 - a 01 ) + cos m -y sin (a Q4 - a Q1 )] + 


^ r IT IMT . f IT , 

+ Qq [-cos m j sm a Q1 + cos sin (m Q y - a^) + 


N>J=! 
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3tt . , Sit . , 3tt ^ .. ■, 

+ cos m sm - ct^ + cos m sm (m o ~ - J 


3a 


rJ IP ==> 0 = q 2 [-cos —■ sin (a^ - + 


02 


+ cos sin (a 03 - a Q2 ) + cos mTr sin (a Q4 - a 0? )] + 


+ Qq [cos m j sin O o j -a Q2 ) + 2 cos ^ sin (m^ - a Q2 ) + 

3tt . , , 3tt . , 3it , 

+ cos m ~J sm (-a Q2 ) + cos m — sin ( m D ~ -a 02 J J 


3 

cia 


03 


IP ==> 0 = q 2 [-cos imr sin (a Q3 - a^) - cos —• sin (a 03 - a Q2 ) + 


+ cos ~ sin Ca Q 4 " a Q3^ + t cos sin ( m 0 7T " a 03^ + 


mu , 


imr . , 3ir , 3 tt . . u _, •> 

+ cos sin (m o —j - a Q;5 ) + cos m — sm (m o j - a 03 J - 


Sir . i 
cos m sm a Q3 J 


^ — IP ==> 0 = q 2 [-cos m ~ sin (a Q4 - a Q1 ) - cos mir sin (a Q4 
u 04 


a Q2 ) 


cos ^ sin (a Q4 - a Q3 )] + Qq [cos sin (m Q — - a Q4 ) + 


m7T 


3tt 


3tt . , Srr tt 

+ cos m — j sin - a 04 ) + cos m sm (m Q ^ " “da-' ' 


04 


7tt . -i 
cos m -j sm a Q4 ] 


For m = 0 and m = 0 , 
o 


0 = q (sin (a Q1 - a Q2 ) + sin (a Q1 - a Q2 ) + sin (a Q1 - a Q2 )} + sin a Q1 
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0 = q {-sin (a 01 - c 02 ] + sin (a Q2 - a^) + sin (a 02 - a^)} + 

+ 4Q sin ctQ 7 

0 = q {-sin (a 01 - a Q3 ) - sin (a Q2 - a Q3 ) + sin (a Q3 - a Q4 ) } + 

+ 4Q sin a^ 3 

0 = a {-sin (a 01 - a Q4 ) - sin (a Q2 - a Q4 ) - sin (a Q3 - « Q4 ) } + 

+ 4Q sin a Q4 

0 = 2q {2 = cos (a Q1 - a Q2 ) + cos (a Q1 - a Q3 ) + cos (a 01 - + 

* COS (o 02 - a 03 ) ♦ cos (a 02 - » 04 ) + cos (a 0J - o 04 )} 

The predicted solution, q = Q and ct^ = = a^ 4 = (2£+1 )tt, 

A = 0, ±1, ±2, . , satisfies the above relations. For this solution, 

th 

the index of performance for the m mode is: 

F 2 

jp _ — 2iL [l - cos [4 + 8 cos my + 4 cos irnr] 

2b 

Hence , ? -? 

F Q 

IP = 0 except when |m| = 4, 12, 20, 28, . . , where IP = 32-^jj — . 

b 

Thus, complete cancellation of the m = 0 mode is expected. Signifi- 
cant sound pressure reduction should be obtained below the cutoff 
frequency for m = 4. Sound pressure level reduction problems can 
be expected at the cutoff frequencies of the radial distributions 

of the I ml =4, 12, 20, , . . modes because the F term in the 
' 1 9 my 

index of performance has singularities at these frequencies. Physically, 
this results in duct resonances of the specified higher order circum- 
ferential modes. Near these frequencies the predicted optimal solution 
no longer holds because a higher order mode propagates, Optimization 
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over the entire set of modes ( I n IP ) would be required to determine 

v m=0 m ^ 

the optimal solution. However, excellent results are predicted for 
f <1913 Hz, the cutoff frequency for the m = 4, y = 0 mode, utiliz- 
ing the predicted solution. This accounts for the good experimental 
results obtained for low frequencies. 

The generation of the m = 4, y = 0 mode is illustrated by the 
circumferential distribution shown in Figure 56, 


2 cancelling waves 
cos m0 for m = 0 



FIGURE 56 

Generation of m = 4 Mode While Cancelling m = 0 Mode 


where t = constant, the symbol • designates the normalized amplitude 
of the source array transducers at time t, and the symbol o desig- 
nates the normalized amplitude of the cancelling array transducers. 
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For m =1 and m = 1. 
o * 

0 = q {-sin (a Q1 - a Q3 )} + /T Q {sin a Q1 - cos a Q1 } 

0 = q {-sin (a Q2 - a Q4 ) > + /T Q {sin a^ 2 - cos a Q9 } 

0 = q { sin (a Q1 - a Q3 ) } + f~2 Q {-sin a Q3 + cos a^} 

0 = q {-sin (a Q2 - a Q4 ) } + /T Q {sin a Q4 + cos a Q4 > 

0 = 2q {2 -cos (a 01 - a Q3 ) - cos (a Q2 - a Q4 ) } + 

+ 2 Q cos a Q1 + sin a - cos a Q2 + sin a Q2 - cos a Q3 

-sin a__ + cos a_. - sin a rt ,} 

03 04 04 

The predicted solution, q = Q and a Q1 = a Q2 = a Q3 - 
3u 

and a Q4 = ~j, satisfies the above relations. For this solution, the 

index of performance for the m^ mode is : 

2 2 
F Q 

IP = ■ {4 (1 - cos imt) + 4 f~2 (cos m - cos m j) } 

b 2 4 4 

Hence, 


IP = 0 except for |m 


IP = 16 


F 2 Q 2 
my x 


3, 5, 11, 13,. 19, 21 


where 


For reference, the IP for q = 0 is nonzero for |m| = 1, 3, 5, 7, . . 
and is given by 



Thus, cancellation difficulties can be expected for frequencies above 
1511 Hz, the cutoff frequency for the m = 3, y = 0 mode. The genera- 
tion of this mode is illustrated by the circumferential distribution 
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shown in Figure 57. 



Figure 57 

Generation of m = 3 Mode While Cancelling the m = 1 Mode 

The symbol • designates the normalized amplitude of the source array 
transducers, and the symbol o designates the normalized amplitude of 
the cancelling array transducers. Examination of Figure 21 , the radial 
distribution of cancellation of the m = 1, y = 0 mode, illustrates 
this phenomenon. The radial distribution for the cancelled condition 
is a high sound pressure level y ~ 0 curve. This can be explained 
by the generation of an m = 3, y = 0 mode above 1511 Hz. 


133 


For 

m 

c 

= 2 

and 

m = 2 

: 






0 = 

2 

q 

{sin 

(“oi 

■ ‘W 

+ sin 

(a oi 

■ 

- sin 

(a 01 

- v 1 

0 = 

2 

q 

{sin 

< a 01 

■ a 02^ 

- sin 

(a 02 

- “03 5 

+ sin 

(a 02 

■ “W* 

H 

O 

2 

q 

{-sin(a 01 

" a 03^ 

+ sin 

' a 02 

- °W 

- sin 

^ a 03 

- “04 )} 

0 = 

2 

q 

{sin 

fa oi 

" °W 

- sin 

( a 02 

- V 

+ sin 

^ a 03 

- “04 )} 

o = 

4q 

{2 - 

COS 

Ca oi - 

a 02> * 

COS 

(“01 - 

i 

< — i 

to 

o 

3 

cos 

(“oi ‘V 


- 

cos 

Ca 02 

" °W 

+ COS 

f “02 

" C W 

- cos 

(“02 



, , . ^ , 37T it 3tt 

The predicted solution, q = Q, and = -y, = y, = -y, 

and cIq^ = y does not satisfy the last equation. The only physically 

reasonable solution is q - 0. Therefore for n = 4 the m = 2 

mode cannot be cancelled. The index of performance of the m^ mode 

for m =2 and q - 0 is : 
o 


IP = 


2 2 
2F Q 
my x 


{(1 + cos mu) - 2 cos 


Hence , 

IP = 0 except for |m| = 2, 6, 10, 14, . . . where 


IP = 8 


2 2 
F Q 
my x 


Because of multiple mode excitation, m Q > 2 will not be considered. 
It is clear that for the generation of these modes, n must be increased 


from n = 4. 
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-The following conclusions can be drawn from this analysis: 

1) There exists a minimum number of transducers, n . , which can 


generate a single mode m^ to a specified frequency f max - This 
number n . can he determined by proper manipulation of the equation 

for the index of performance developed in this section. 

2) For a cancelling array interacting with a source array, as in 

our case, there is an additional restriction for a cancellable condition. 

There exists a minimum number of transducers, n . , which can produce 

mm 

a cancelled condition for the m Q mode to a specified frequency f max 

without the generation of higher order modes. 

3) For the case of two interacting arrays, the second array can 
only impose further restrictions on the higher order modes which can be 

i 

generated. Therefore, n . < n . . 

° min — min 

The number n . can be determined from the index of performance 
min 


for q = 0. 


IP 


F 2 

my 


n n 


2b k 


2 2 cos m (6 ok - 0 ) cos (* k - Q"> 

= 1 j = l J J 


F 2 Q 2 n-1 - 2Trm 

= — {n + 2 2 (n - j) cos — — - j cos — — j} (3.3-C) 

2b ^ j=l 


Hence, n is increased from a trial value until IP = 0 for all m 
except m Q whose y = 0 cutoff frequency is below f max - The final 

value of n is n . . 

min 
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The generation and cancelling of submodes of m Q by a single array 

must also be considered. A large portion of acoustic power can be lost 
by this phenomenon. For example, consider the case of n = 4, m Q = 2 . 

As illustrated below, acoustic energy can be wasted in the generation 
of two cancelling plane waves as shown in Figure 58, where 


cancelling plane waves 



= 2 
0 


FIGURE 58 

Generation of Cancelling Submodes 


t = constant and the symbol • designates the normalized amplitude of 
the transducers at time t. Further testing would be required to verify 
the existence of this phenomenon. This particular case illustrates one 
additional problem. The relative positioning and phasing of trans- 
ducers makes impossible the cancellation by an array displaced circum- 
ferentially — radians. The cancelling array only provides a bias 
for the direction of spin and increases the sound pressure level. So, 
the additional restriction n > 2m Q is required. 
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The range of m, y, and frequency over which effective cancella- 
tion can be obtained can now be cataloged. Table 8 shows the 
results for n = 4 , Upper frequency limits have been set approximately 
100 Hz below the cutoff frequency of the higher order mode subject to 
excitation. This margin assures that the decay rate in the axial 
direction of the higher order mode is great enough to have little 
effect on the magnitude of sound pressure level reduction. 

TABLE 8 

m, y , and f Range for Cancellation with n = 4 


F range 


O 

il 

3. 

y = i 

y = 2 

m = 0 
0 

0 

to 1378 

1378 to 1800 

- 

m o = 1 

662 

to 1400 

- 

- 

m o = 2 


- 

- 

- 


Examination of the experimental results of Chapter II substantiates the 
above figures. The actual results are slightly higher as a result of 
the duct resonant frequencies being slightly higher than predicted by 
theory . 

The analysis of the wave equation has shown that the frequency 
range of propagation of a single, specified circumferential mode gen- 
erated by a single array is finite and dependent upon the number of 
transducers. A circular transducer array given by equations 1.3-1 
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through 1.3-4 where n is finite and m is specified (m o ) can 

excite an infinite set of circumferential modes, {Mi. The magnitudes 
of the elements of {M} are given by the |m| which yield a nonzero 
solution to the sound pressure level index of performance for a single 
array: ~ - 

F“ Q n-1 0mT , 2irm 

IP = JUiL — {n + 2 £ (n - j) cos ^ j cos j) (3.3-C) 

2b / j=l 

For a given m o , the spacing of the elements of {M} generally increases 
as n is increased. This is illustrated in Table 9 for m Q = 1. 


TABLE 9 

{ |M(n) | } for m Q = 1 


m » 1 

o ! 

{ )m| } 

n - 1 

n > 1 for equation 3.3-B to hold 

n = 2 

{1, 3, 5, 7, 9, . . .} 

n = 3 

{1, 2, 4, 5, 7, 8, 10, 11, . . .} 

n = 4 

{1, 3, 5, 7, 9, . . .} 

n = 5 

{1, 4, 6, 9, 11, 14, 16, ... } 

n = 6 

{1, 5, 7, 11, 13, 17, 19, . . .} 

n = 7 

(1, 6, 8, 13, 15, 20, 22, 27, 29, . 

00 

11 

{1, 7, 9, 15, 17, 23, 25, . . .} 


The information in Table 9 was obtained by evaluation of the bracketed 
term in equation 3.3-B. A simpler approach yielding more information is 
available. Tyler and Sofrin (4) have derived a simple equation which 
yields the set of circumferential mode numbers which can be excited by 
rotor-stator interraction in a fan. The value of m is restricted to: 

I 

m = n B + kv 
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where 

m = circumferential mode number 
n'= harmonic number 
B = number of blades 
k = . . • "1, 0, 1, . • • 
v = number of vanes 

A similar equation can be derived by analogy for a transducer array. 

t 

The array is excited only at frequency w, so the harmonic number n 
is restricted to the fundamental. The transducers are stationary 
sources, so v is set equal to the number of transducers, n. The 
value of B is selected so that the individual rotor-stator inter- 
actions of the blade-vane system coincide with the positive maximum 
of the sinusoidal transducer outputs. Therefore B is simply m Q . 
The rotor of the fan would have the same angular frequency as the 
acoustic frequency of the transducers. The equation can now be re- 
written: 

. (3.3-E) 

m = m + to 
o 

To illustrate the effectiveness of this equation, consider the 
case of m Q = 1. Equation 3.3-E yields the sign of m (recall that 
a negative sign denotes spin in the opposite direction of the rotor) 
So the set {M> can be completely specified. Table 10 contains 


the new information. 
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TABLE 10 

(M(n)} for m o = 1 


m = 1 
0 

{M} 

rH 

il 

{0, + 1 , +2 , +3, +4, . « .} 

3 

11 

{+1, +3, +5, +7, . . .} 

n = 3 

{T, -2, 4, -5, 7, -8, . . .} 

n = 4 

{1, -3, 5, -7, 9, -11, . . .) 

n 55 5 

{1, -4, 6, -9, 11, -14, . . .} 

n - 6 

(1, -5, 7, -11, 13, -17, . . .} 

n = 7 

{1, -6, 8, -13, 15, -20,. . . .} 

CO 

ii 

£ 

{1, -7, 9, -15, 17, -23, . . .} 


Equation 3.3-E provides a very simple means for obtaining the minimum 
number of transducers in a single array necessary to generate a single, 
specified circumferential mode over a specified frequency range. 

The blade-vane analogy can also be applied to the two array system 

with the cancelling array set for maximum cancellation by equations 

1.3-7 and 1.3-8. The analysis of the transducer-duct system leads to 

a model of a rotor-stator source with 2n vanes, (n + m ) blades, and 

an angular rotor frequency of — — - — times the acoustic frequency of 

o . 

the transducers. The empirical relation yielding the circumferential 
modes which form the set {M} is then given by: 

m = m Q + (2k + 1) n (3.3-F) 

The sets of {M} for the four transducer array set for cancellation can 
now be easily specified: 
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TABLE 11 

{Mfm o )} for the 4 Transducer Arrays Set for Cancellation 



{M} 

3 ! 

O 

11 

o 

{+4, +12, +20, +28, . . .} 

m ® 1 

{-3, 5, -11, 13, -19, 21, . . .} 

0 


m - 2 

{-2, 6, -10, 14, -18, 22, . . .} 

0 


m =3 

{-1, 7, -9, 15, -17, 23, . . .} 

0 



The analysis of the sound field is hence complete. To complete the 
mathematical model of the system to be controlled, a measurement system 
must be produced which can generate a sound pressure level index of 
performance. Also, the transfer function for the transducers must be 
obtained. 

An ideal measurement system for determining an overall sound pres- 
sure level for a sound field of unknown distribution would yield a 
function of the form: 

SPL(t) = ^ f r / Q f z SPL (r,0, z,t) dzd6dr 

Trr z 

which is similar to the index of performance which has been used. With 
a finite number of stationary probes, a sound pressure level indicator 
could be defined by: 

1 w 

SPL (t) = i i Z 1 e i SPL.(t) 

where w = number of probes 

6^ = weighting factor related to position of probe 
SPL^ = SPL at i position 



141 


= 20 log Pi -20 log P ref where Pi is the RMS value of Pi 
determined by equation (3.3-A) 

Such a system is physically relizable and can be designed to be 
relatively insensitive to the point cancellation effect. This form 
shall be assumed for the development of the controller. 

The transducers must also be modeled. A ducted transducer has 
the following properties: 

1) SPL.. = ki Cf) + k 2 (f) (20 log V. -20 log V ref) 

where k^f) and k 2 (f) are frequency-dependent parameters 

v. = electrical voltage across transducer. 

1 

"th 

SPL. = SPL near the j transducer. 

3 

2) a . = 9. + k_(f) 

J 03 3 3 

where a . = acoustic phase 

oj 

0. = electrical phase 

3 

kg = frequency dependent parameter 
The source strength q can be related to SPL^ : 
q = f (SPL.). 

The plant can be modeled as shown in Figure 59. The diffi- 
culties associated with this plant can be easily seen by this model. 

As found in Chapter II, the functional relationships of k^k^ and 
k^ to frequency are too complex to be described analytically and 
provide an unreasonable function generation problem. In reality, 

k (f), k 7 (f) and k,(f) also vary with the transducer and its location, 
12 d 

though the frequency dependence is quite similar from point to point. 



With the simulated source, Q and cj> o j are also dependent upon a 
similar set of k's. In addition the relation for the pressure field 
is extremely complex. 




? f (kj (f) + k 2 Cf) (20 log v^ -20 log v yef ) 


a . = 0. + k_(f) 
oi 1 3^ 


spinning wave 


spinning wave 


CO co 


R[P(r,0,z,t) ] = E E e J (k r ) J (k r) • 

1 J m=0 y=0 my m m v my o J m v my J 

n 

* E (Q cos (d> . + y z - cot) cos (m0 - m0 ,) 
x vr oj 'my oj' 


+ q cos (a . + y ,z - cot) cos (m0 - m0 . - — -) } 
n oj 'mil o j n 


combined acoustic field 



Mathematical Plant Model 
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Conventional optimal control techniques are not appropriate because 
of the complexity of the plant. Also, the sound field equations have 
been developed for a specified source, a transducer array. So, a 
conventional controller would be restricted to this type of simulated 
fan source. The application of adaptive control techniques would thus 
be appropriate and desirable. Then the SMS cancelling array could be 
applied to a variety of acoustic noise sources without modification. 

The use of adaptive control had been proposed in the original 
research proposal. Because of the complexity of the transducer trans- 
fer function, the decision to design an adaptive controller was not 
made until the completion of experimentation. At this point, suf- 
ficient evidence had been gathered to support the adaptive design 
approach. The analysis presented in this section was developed 
parallel to the development of the controller. In addition to pro- 
viding an understanding of the plant and providing background for 
future development, it gives further support to the original design 
proposal . 
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3.4 Synthesis of Adaptive Control Scheme 

The model presented in the previous section does not entirely 
encompass the behavior of the plant. The relationship for P(r,0,z,t) 
assumes symmetric waveforms due to a symmetric duct and ideal trans- 
ducer outputs. Thus for the predicted solution and a Q must be of 
the form: 

3 .(t) = q(t) [1,1,1,1] T (3.4-D 

a Q (t) = a D ( t ) + [0, m j, mir, m^-] T (3.4-2) 

where a is in radians, 
o 

This also implies symmetric acoustic loading of the transducers so that 
k^(f), k^f), and k^(f) are each uniformly alike for all transducers. 

These assumptions reduce the dimensionality of the system from 
eight to two. So, the design of a controller based on these assumptions 
is an attractive proposition. Its justification is twofold: 

1) the scope of concepts which can be effectively investigated 
would not be seriously limited by the suboptimal manipulation 
of the transducer array. 

2) a two dimensional controller can serve as a rapid convergence 
initial stage of a multi-stage controller for optimal cancel- 
lation. 

The two dimensional system offers obvious design, construction, and 
adaption advantages. So, the decision to design a 2-dimensional 
adaptive controller was made. 

The plant controlled variables for the ideal plant can be defined 
by: 
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x T (t) « [x^t), x 2 (t)] = [v(t) ,e(t) ] 

such that vj(t) = v(t) [ 1,1,1 ,1] 
and 0 T (t) = [0, j, tt, + [9, 0,e,e] 
where v(t) = electrical voltage and 6(t) = electrical phase. The 
measurement system formulated in section 3.3, 


1 w 

SPL(t) = “ £ &. SPL (t) 

i=l 

where w = number of probes 

fL = weighting factor related to position of probe 

+■ Vi 

SPL^(t) = SPL at i position, 

provides a meaningful index of performance or cost function. It shall 
be designated as J(x(t),t), indicating its dependence on the plant 
controlled variables and time. Maximum cancellation occurs when 
J(x(t),t) is minimized. The optimum inputs providing J min shall be 
designated by X*(t) = [X 1 *(t), X 2 *(t)]. 

The cost function with an arbitrarily located minimum would appear 
as in Figure 58. This applies for the predicted form of the solution 
when the system is operating under modal conditions satisfied by n n) ^ T1 
(see section 3.3). The following features are important: 

1) 2-dimensional plant 

2) unimodal cost function surface 

3) time varying minimum 

4) probably a noisy surface because of the SPL measurement 


algorithm. 



Figure 60 

A Hypothetical Realization of 
the Cost Function Surface 






X|(t) = V(t) 

UNITS = 20 LOG ( V 1 .5 ) 
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An adaptive controller would use a search technique to locate the 
minimum and follow its time variation. There are many techniques 
applicable to our particular problem. Those considered for application 
are presented in the next chapter. 
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IV Adaptive Optimization Procedures for Multi-Dimensional, Unimodal, 
Noisy Systems 

A literature search yielded several methods satisfying the multi- 
dimensional, unimodal, noisy signal requirements. McMurtry (11) has 
cataloged a representative set of these and other adaptive methods. 
Because of his organized approach to the presentation of method 
algorithms, his work serves as a basis for this chapter. 


4.1 Method of Steepest Descent 

Methods of steepest descent converge to the minimum of a surface 
by moving in the direction of steepest slope. A deterministic surface, 
J(x), is assumed. But, adequate performance can often be obtained on 
systems contaiminated by noise. 

The method is discrete, and determines the gradient of the function 
J (x) at each trial. The next step is made in the negative gradient 
direction. Given the gradient vector J (_x) as 


V 

X 


J(x) = [ 




T 

9J_1 * 

3x J 


n 


(4.1-1) 


the algorithm is 


x, .. = x. - c. V J (x) (4.1-2) 

— k + 1 -k k x — 

where c^_ can be a constant or variable gain sequence. An adaptive 
adjustment procedure for c^ has been proposed by Eveleigh. Increased 
rate of convergence and decreased steady state error at x* (where 
J(x*) = 0) can be expected from such a scheme. 



ISO 


A continuous controller based on this method could utilize the 
algorithm: 

9J 


dx. 

1 


dt 


= - K 


9x 


i = 1,2,... ,n. 


(4.1-3) 


The gradient can be determined continuously by a perturbation technique. 


4.2 Stochastic Approximation 

Stochastic approximation methods, the stochastic analogs of the 
steepest descent method, may be employed to optimize systems corrupted 
by noise. This method can be applied to time- invariant systems. To 
establish notation, the problem is stated as follows: it is desired 

to minimize an unknown, deterministic function J (x) where x = [x^,...,x ]' 
are controlled variables. Contaminating the system output is a zero 
mean random process N(x). Thus E[J(x) + N(x)] = J (x) . 

Consider the minimization of J (x) by finding the zero of V J (x) 

X 

in the presence of noise. Let Y(x) be the noisy value of the gradient 
function where E[Y(x)] = J (x) . Thus Y_(x) is a random vector. The 
algorithm is: 

2k +1 + 2k - fi k 

where 6^ is a variable gain sequence. It can be shown that when the 
following conditions are satisfied 


1) \ > o 


2 ) Z &. = « 
k-i k 


Z < <*> 

K=1 k 


(4.2-2) 


3 ) 



5) Vx_ J (x) behaves linearly for x^ near x* 


the algorithm converges to x* (where V x J(x*) = 0) in mean square 
and with probability one. 

Kiefer - Wolfowitz methods enable the determination of the extremum 
of J (x ) when Y^x) is not available. For M(x) = J (x) + N(x), a gradient 
can be determined by 


n 

V y M(x) = Z e. [M(x^ + c k e^) - Mf^ - c R e..)] 


j=l j 2c k 


(4.2-3) 


. th 


where e^ are unit vectors in the j direction and c k is the distance 
on either side of Xj, at which measurements are made. The algorithm 
is : 


’ 2k - \ M(x) 

It can be shown for the Dvoretsky conditions: 
lim 

1) k^» = 0 

lim 

2) k-*» c k = 0 
lim n 

3) n-* 30 Z 5, = 00 

k=l 

lim n 6, 

4) n-*» 2 ( — -) 2 < °° 

k=l C k 


(4.2-4) 


5 ) l J(; ik + l^ ' J( ^k^l < A l^k+i - 2^1 + B < 00 
that the algorithm converges to x* (where (x*) = 0) in mean square 
and with probability one* Convergence acceleration techniques are 


available. 
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4.3 Adaptive Step Size Random Search 

Schumer and Steiglitz developed a random search method that 

uses the optimum step size at each step. The proposed adaptive step 

size random search has been shown to approximate the optimum step size 

random search (OSSRS) for a hyperspherical surface only. 

The algorithm is as follows: two steps are made in random di- 

th 

rections of length s^ and s^ (1+a), 1 > a > 0, on the k iteration. 

Sj^ is set equal to the step size yielding the greatest improvement. 

If there is no improvement after a specified number of iterations, the 

step size is reduced. To ensure that s, does not become too small, after 

! k 

a large number of iterations a is made much greater than one. This 
enables comparing s^ with a very large step. 

4.4 Matyas' Method 

Matyas considered adaptive optimization by the random variation of 
parameter magnitude and direction. Convergence of a random sequence of 
parameters for a stationary unimodal surface has been proven. The con- 
cept of directional adaptation, requiring memory, provides an algorithm 
for an increased rate of convergence: 

4*i ' 4 * 4.1 * T k.i 4.1 = 4 * 4.i 

where T, - is a variable correlation rtiatrix, C, , is a random vector 
k+1 Hc+l 

selected from a multivariate normal distribution with zero mean and 
unit correlation matrix, ^ +1 = d^ +1 + T R+1 C k+1> and d^ is given by: 


d. , = [c (x.) + c. (x. ) 1 
-k+1 o -k^ lMc 1 


k-1 

d, + Z h. , .6. 
-k l ,k- l—i 


+ c i ( 4> \ 4 


(4.4-2) 
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where 

0 1 C 0 Q^) 1 X > > °> and + c l^) > 1 if s k = 0 

0 1 Cq^) 1 1 > 1 °> and ICqC^) + Cxj.) j <i if <$ k = 1 

(4.4-3) 

and 

k-1 

h. . . « H 6.. 

i.k-1 j,! J 

Although requiring greater computing capability than the simpler methods 
already mentioned, the method can be easily implemented on a digital 
computer. 

4.5 Gradient Biased Random Search 

Pensa used the noisy gradient measurement (equation 4.2-3) as 
the bias vector in equation 4.4-1 of Matyas' method. The rate 

of convergence on a noisy unimodal surface is optimized by specifying 
the covariance matrix . Convergence has been proven for a noisy 


unimodal surface. 
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V Perturbation Types Adaptive Controller 
5.1 Perturbation Signals 

The gradient of a cost function surface can be determined con- 
tinuously by the perturbation, or small amplitude oscillation about 
the current value, of the input signals to the system. This in effect 
increases the dimension of the output by producing information about 
the system as it responds to the perturbations. Thus, for a single 
output, n-dimensional system with n perturbated inputs, the output 
can contain up to n+1 simultaneous, separable signals. Filtering 
provides the signal separation. 

With signal perturbation, a continuous -time gradient technique can 
be applied. This would require only a relatively simple continuous 
controller, since discrete operations, memory, and logic, requiring a 
digital computer type controller, would not be needed. With conditions 
on the perturbation amplitudes to be discussed later, this technique 
can work well on noisy cost function surfaces. Therefore, it will 
be investigated in depth. 

The plant to be controlled is a high frequency electrical-acousti- 
cal system. For practical experimental purposes, the acoustic frequency 
can be assumed to be greater than 600 Hz. So, equilibrium conditions 
and good signal separation can be obtained for perturbations to the 
system of less than about 100 Hz. Voltage perturbation will be realized 
as amplitude modulation of the high frequency signal. Phase perturbation 
will simply be modulation of the signal phase controller. 

The following obtained from a perturbation signal can be analyzed 


as follows: 
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Let J(X. ) be one dimension of a smooth, unknown cost function where 
X is slowly time varying for the time interval of investigation. 


jc>y 



FIGURE 61 

One Dimensional Perturbation 


The perturbed signal is: = x iopt + + A sin io 1 t 


Let x k = X i - X opf Then \ = x m + A Sin “l 1 


(5.1-1) 

( 5 . 1 - 2 ) 


Use a polynomial approximation for the unknown cost function in the 
region x-c<x^<x+ c, c>A: 


J * a 0 * Vic + a 2 x k * 


The slope of J at is: 


n 

+ a x. 
n k 


( 5 * 1 - 3 ] 


3J 

3x, 


x i = 

k m 


= a n + 2a 0 x + 3a + . . . + 
1 2 m o m 


n m 


(5.1-4) 
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The output J (t) due to the input + A sin t is: 

2 

J(t) = a A + a, (x + A sin go, t) + a„ (x + A sin go, t) 

v ' 0 1 ^ m 1 ' 2 v m 1 

(5.1-5) 

+ a„ (x + A sin go, t)^ + . . . + a (x + A sin to, t) n 

3 m 1 n m 1 

2 

' = a- + a, (x + A sin w. t) + a„ (x + 2A sin to, t + 

01m 1 2 v m 1 

12 3 2 

+ ■=■ A (1 - cos 2 go, t) ) + a, (x + 3x A sin go, t + 

2 13mm J 

+ y xA^ (1 - cos 2 go^ t) + A? (3 sin GO^ t - sin 3 to^ t) ) 

4 3 2 2 

+ a, (x + 4x A sin go, t + 3x A (1 - cos to, t) + 

4 v m m 1 1 

3 14 

+ x A (3 sin to, t - sin 3 go, t) + - 5 - A (3-4 cos 2 to, t + 

m ^ 1 18 1 

+ cos 4 to, t)) + . . . + a (x + A sin to, t) n 
1 n m 1 


The component of J (t) which varies at frequency u^, which can be obtained 
by ideal filtering, is: 

2 3 2 3 2 

y = (a, + 2a„ x + 3a„ x + -r a_ A + 4a. x + a. x A +,..)• 

7 ^ 1 2m 3m 43 . 4m 4m 

* A sin GOj t (5. 1-6) 


The amplitude of y is: 


A 

y 



+ . . .) 

(5.1-7) 
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Defining the parameter S to be 


S - (»! * 2=2 \ * 3a 3 x m * T a 3 A * * 4a 4 % * a 4 %, * ' ' 0 ' 


it can be seen by comparison with Equation 5.1-4 that the slope of J 

at x, = x can be approximated by S: 
km * 1 


3J 

^k' 


x 


m 


S 


(5.1-8) 


The error is given by: 

e = 4 a, + a. x A^ + . . . (5.1-9) 

3 3 4 m 

The error is due to the higher order terms of the polynomial approxima- 
tion and is negligible for » A. As x^ -*■ 0 however, A becomes 

a contributing factor to the steady state error. 

For m independent perturbed inputs (uk jt for all j,i and 

£=1,2,3. . .), m components of the cost function gradient can be 
determined. The response of the system will be limited by the time 
required to determine A^.. Error is accumulative only over the interval 

during which A^ is determined. 

To minimize error in the determination of the gradient due to 
higher order components of the cost function surface and to minimize 
steady state error, it is desirable to minimize the perturbation ampli- 
tude, A. But, system irregularities, nonlinearities, and noise 
constrain the minimum value of A. It is therefore important to deter- 
mine an A . as a function of the parameters of the non-ideal system 
mm 


characteristics . 



158 


Output signal separability is essential for complete perturbation 
measurability. So, each filtered perturbation output signal, y t > must 
be independent of all other y., j * i. Thus, a 1-dimensi.onal analysis 
can be applied to each of the perturbed system dimensions. Shown below 
is a realization of J versus x i with random noise present for the 
i^ 1 dimension at f = f^ and x^, j j- i, fixed. 


J 



FIGURE 62 


til 

Realization of Cost Function in i Dimension 

Of importance is the departure from a smooth, ideal function. This super- 
imposed random process which is dependent upon x^ shall be denoted by 
the symbol N (X.). This random process was detected experimentally, 
and is dependent upon the acoustic frequency and the values of the 
other input parameters. The random process arises due to the point 
cancellation effect on a SPL measurement network. Individual micro- 
phone probes are effected independently because of localized distor- 
tions of the sound field due to an out -of ground duct, variation of 
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transducer- to-transducer response characteristics, etc. The measure- 
ment network combines these signals to obtain the "noisy" performance 
indicator, J. The output is dependent upon the acoustic frequency, 
which is a continuous variable from 600 Hz. to 6000 Hz. , and the other 

parameter, X., which is continuous from 0 to X. . Thus, there 
y j max 

exists an infinite number of realizations, making deterministic 
analysis inapplicable. 

We have for a single dimension 



JCX^ 

= f(X.) + N(X.) 

where 

X. = 
i 

x. + A sin to t 
10 


£(x i ) 

= smooth deterministic function 


N(X.) 

= random component 


Since we are dealing with only 1 dimension, the i subscript will be 
temporarily dropped. The following qualifications will be used to 
establish a mathematical model for N(x). 

1) The rate of change of the acoustic frequency is slow 
enough that the realization of N(x) remains unchanged for 
the time period for which J(x) is being examined. This 
makes n(X) slowly time varying in control systems context. 

2) Each acoustic frequency has a unique realization of N(x) and 
the selection of the acoustic frequency is not a deterministic 
process. 

N(X^) is independent of N(Xj) for all i, j and i ^ j . 


3) 
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It is therefore appropriate to model N(x) as a random process {N} 
such that 

1) P N (X) - j-ft e 2 °N (Gaussian) 

where c?, T is the standard deviation of the random process 
N 

2) (N(X)} is white and of zero mean. 

Thus : 

E [J (X) ] » E[f(X)3 + E [N (X) ] 

= f(X) (5.1-11) 

System nonlinearities may cause a vertical shift of the transfer 
(cost) function. Let the magnitude of this shift be Define: 

g(X) = f(X) + l. (5.1-12) 

g(X) can be used as a new cost function surface without effecting the 
position of the minimum or the operation of the controller. Only the 
dc component associated with X q is changed, which is filtered out. 
Now, 

J(X) = g(X) + N (X) (5.1-13) 

Thus , 

E[J(X)] = g(X) (5.1-14) 

Therefore, the expected value of an output signal is that obtained for 
a smooth transfer (or cost) function plus a constant. The analysis 
for a smooth, or ideal, cost function then applies to the expected 


signal . 
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Information is obtained from the <h. frequency component of the 
output, where on is the perturbation frequency in the i 1 * 1 dimension. 

The processing of this signal introduces a bias due to N(X). In 
particular, say the RMS value of the signal is used to determine its 
amplitude. 


E[J 2 (X)] = E[g(X) + N (X) ] 2 

= E [g 2 (X) + 2g (X)N(X) + N 2 (X)] 
= g 2 (X) + 0 + a N 2 


From (SI. -6) and (5.1-8), the filtered output from a smooth cost func- 
tion surface at the perturbation frequency is: 


dJ 


y ■ A ax sin V - J'i 


dJ 


Including the noise, y^ = A ^ sin to^t + N(X, co^t) 


(5. l-8a) 
(5.1-16) 


where N(X,Ujt) = N(X(t) + A sin to t) 


w = (jo. 


E[yi 2 (t)] = E[A 2 (g) 2 sin 2 u^t + 2A sin to 1 1 N(X,cot) + N 2 (X,Lo.,t)] 


dX 


. 2 ,dJ, 2 .2 „ A 2, , 

= A (^) sin Ujt + 0 + a N (Wj) 


(5.1-17) 
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The mean value can be determined: 

y l 2(t) = ? ; o y l 2ct) dt 

X 7T 

it [A 2 (— ) 2 sin 2 o^t + a ^ 2 (u) ] dWjt 

= “ {A 2 (-~) 2 Ujt - j sin 2 u 1 t) + (uijt) a N 2 (a) 1 )] 71 

o 

= ^HX + °N (5.1-18) 



or A > 
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can be determined experimentally. So for N e specified, 


A . = 


10 OjjCwj) 


min 3J 


(5.1-21) 


ax. /n 

i e 


3 J 

For a constant amplitude perturbation system, the value of f° r which 

9 J 

\iin v ^°^ ates ^ is specified. Thus, letting S = establishes 


A . . 
min 

The a./w) obtained as a result of the "variation" of the cost 
N 

function surface is a parameter of noise at frequency The output 

signal will also be contaminated by signal noise in the usual sense. 
There exists a random process, presumably Gaussian, of zero mean and 

I 

standard deviation cr^ (la^) which contains all signal noise. The mini- 
mum perturbation amplitude is then specified by: 


A . 


mm 


10 a N (ojj) 
S /U 

e 


( 5 . 1 - 22 ) 


5.2 Perturbation Circle 

For a 2-dimensional perturbated input, single output system, 
filtering of the input yields the following information: 


3J 

ax. 

i 


(t) 


V ct) 
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where and are the amplitudes of the input perturbations 

and Ay^Ct) and A^Ct) are the respective outputs. 

To simplify the control system, the concept of a normalized 
system is introduced. Here, all quantities in the controller are 
proportional to those in the actual system by the proportionality 
constants a and g. But the proportionality constants are selected 
so that corresponding quantities in the controller have the same magni- 
tude. For example, for the control variables and X^, let 

X L max (controller) = aX^ max (Actual) = K volts 

X2 max (controller) = g X ? max (actual) = K volts 

Thus the controller inputs and outputs would be amplified or reduced, 
depending on the constants a and g. 

Using the proportionality constants a and g, we can let 
Aj = A^ = A. Then, 



The angle of the negative gradients direction of J on the (Xj,X 0 ) 

* 

plane with respect to the positive axis, y , can be obtained from 
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the definition of the gradient: 


3J 

ax. 


sin y = 


,aj ,2 ,aj ,2 

<5xp> 


3J 

ax. 


cos y = 


,3J 2 ,3J 2 

l 3X 1 J L 9X 2 J 


- V (t) 


V 2ft) * A y2 2 "> 


yi 


A y2 2(t > + A yl 2(t) 


(5.2-la) 


(5.2-lb) 


For a gradient method, the direction of the negative gradient, y , 

is the desired direction of movement to converge on the minimum. 

* * 

If sin y and cos y ate to be found continuously as a function 
of time, the system must be perturbed simultaneously in the and 


directions. For simplicity, it would be desirable to develop a system 
utilizing a single perturbation frequency. In addition, perturbation of 
both inputs at a single maximum frequency will reduce the system time 
response. This could be accomplished by a perturbation "circle" which 
could be created by the inputs : 


X 1 (t) = Xj(t) + A cos 0)jt (5. 2-2a) 

X 2 (t) = x 2 (t) + A sin uijt (5. 2-2b) 

System time delays can be modeled by: 

*l(t) = Xj(t) + A cos (to 1 1 -B 1 ) (5. 2-3a) 

X 2 (t) = X 7 (t) + A sin (uijt -3 2 ) 


(5. 2 -3b) 
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As | $ - BJ the perturbation circle collapses to a single dimension 

perturbation, as illustrated in Figure 63. 



FIGURE 63 

Two Dimensional, Single Frequency Perturbation 

This would result in incomplete observability, and thus, incomplete 
convergence. The phase angles and £ 2 can be measured. With 

calibrated time delays on the perturbation signals, ^ and B 2 

can be increased to 2mr, n = 1, 2, 3, . . ., which has the effect 
of eliminating the time delays from the perturbation analysis. 
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For analysis purposes, assume the X^, X2 coordinate system has 

been centered on the cost function minimum. Then, in the region of 
(x^, x^) , the following polynomial approximation applies: 


J = 


2 

Z 

s=l 


sO 


+ a 


si 


X + 

s 


2 • 
a ^ x + 
s2 s 


+ a x 
sn s 


(5.2-4) 


The output J(t) due to the input 
X(t) = x^ + A cos ujjt 

= x 2 + A sin co^t 


is : 


9 

J(t) = a 1{ j + a ^ (x 1 + A cos Wjt) + a 1 ^(x 1 + A cos Wjt) + 

3 

+ a 13 (Xj^ + A cos u^tj + . . . + a ln (x 1 + A cos Wjt) 11 + 

2 

+ a 2Q + a 21 (x 2 + A sin Wjt) + a 22 (x 2 + A sin Wjt) + 

3 Tl 

+ a 23 (x 2 + A sin u^t) + . . . + a 2n (x 2 + A sin o^t) 

= (a 1Q + a 2Q ) + a^ (Xj + A cos o^t) + a 21 (x 9 + A sin ojjt) + 

2 12 

+ a 12 [x 1 + 2 A cos w t + y A (1 + cos 2 t)] + 

2 12 

+ a 22 [x 2 + 2 A x 2 sin w^t + -j A (1 - cos 2 w^t)] + 

3 2 3 2 

+ a^ 3 (x^ + 3 x^ A cos ai^t + - A (1 + cos 2 ai^t) + 

13 3 2 

+ jA (3 cos to^t + cos 3 Wjt)) + a 23 |x ? + 3 x 9 A sin w-t + 
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3 2 13 

+ j A ( 1 - cos 2 a) t) + A (3 sin to^t - sin 3 w^t)] + 

4 3 9 2 

+ a^ 4 [x x + 4 x^ A cos to^t + 3 x^“ A (1 + cos 2 C 0 jt) + 

3 14 

+ Xj A (3 cos to^t + cos 3 tGjt) + (3 + 4 cos 2 w^t + 

4 3' 

+ cos 4 to^t)] + C x 2 + 4 x 2 A sin uijt + 

2 2 3 

+ 3 x^ A (1 - cos 2 o»jt) + x 2 A (3 sin co^t - sin 3 co^t) + 

14 

+ -g- A (3-4 cos 2 to^t + cos 4 to^t)] + . . . + 

+ a ln ( x 1 + A cos Wjt) 11 + a 2n (x 2 + A sin u)jt) n (5.2-5) 

The component o£ J(t) at frequency , which can be obtained by 
filtering, is: 

y(t) = (a n + 2a 12 x 1 + 3a 13 x ^ + ~ a 13 A 2 + 4a 14 Xj. A 2 + . . .) • 

* A cos u> lt + (a 21 + 2a 22 x 2 + 3a 23 x 2 2 + | &23 A 2 + 4a 24 x 2 3 + 


+ a 24 x 2 A + . . .) A sin oo^t = (A^ cos to^t + A^ 2 sin co^t) 


(5.2-6) 


A ^ and A^ 2 are proportional to the slopes in the X^ and X^, direc- 
tion, as shown in Section 5.1. 

The single separation and analysis is illustrated by the following 
controller operation outline: 



The controlled variables are given by: 


X(t) = 


x^t) + A cos a^t 
x 2 (t) + A sin ajjt 


Let X(t) = k J"** where k is the amplitude of X(t) due to the sinusoidal 
perturbation. In vector form, X(t) is given by 



* 

k cos y 


•k 

k sin y 


The ideal filtered output is: 

y (t) = A y (t) sin (o^t + <f>) 

= A yl (t) cos u>t + A y2 (t) sin o^t 
Multiply Y(t) by 1* cos to^t to obtain: 

u x Ct) = y(t) • cos (jJjt = jA yl (t) (1 + cos 2 w 1 t) + 

+ j A y2 (t) sin 2 w l t 

Filtering u^Ct) with a low pass filter yields: 

z l Ct) ' l A yl (t) 

Multiplying y(t) by 1 * sin a^t yields: 

u 2 Ct) = jA yl Ct) sin 2 cojt + j A y 2 Ct) C l - cos 2 Wjt) 

Filtering u 2 (t) with a low pass filter yields: 

1 
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Equations (5.2-la) and 5.2-lb) are 


sin y (t) = 


-Vt) 


yy * A,%ct) 


* 2 r 

y2 


cos Y (t) * 






2 2 

A (t) = A . (t) + A _(t) , which is the amplitude of y(t) . The full 

y yi yz 

wave rectified, low pass filtered signal y(t) yields 


z 3 (t) = y A y (t) 


Then , 

* 

sin y (t) 


z 7 (t) 

W 


* 

cos y (f) 


zyt) 

^Tt) 


So, 


X = 


K 


Zj(t) 

z 7 (t) 

z 3 (t) 


and 


X(t) = X (t) dt 
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To study the sensitivity of the system to small variations in time 
delays, let 


X(t) 


x^(t) + A cos (u^t - a^) 
x 2 (t) + A sin (a^t - a 2 ) 


where and a ? are small (less than .1 radians). The filtered 

ouput of the system at frequency due to X_(t) would be: 


y(t) = A^ cos (tojt - a 1 ) + A y2 sin (a^t - a 2 ) 


y2 


= A yl [cos a^t cos a 1 + sin Wjt sin c^J + 
+ Ay 2 [ sin a)^t cos a ? - cos u^t sin a 2 ] 
- [cos a 1 A , - sin a 2 A y2 ] cos u^t + 


+ [sin A j + cos a 2 A y2 ) sin u^t 


y2 


S l A yl “ a 2 A y 2^ C0S “l* + t a l A yl + A y2 ] sin “l* 

The two coefficients obtained by the signal separation scheme would be: 


A .=* A , -ou A „ 
yl yl 2 y2 

t 

A = A , + a, A . 
y2 y2 1 yl 

* * 

The control system would generate sin y and y directly from these 
coefficients. For this error analysis, it will be useful to determine 
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| -y f - y*\, the change in the indicated negative gradient direction due 
to and a 2 . Y 1 is the aberrant value of the negative gradient 

direction. For illustration purposes, and a 2 will be examined 

separately. 


y2 

tan y* = T* — 

yi 

A' 2 

tan Y' = 

yi 


A v2 : Vyl 

A yl - a 2 A y2 


tan Y 1 - 


A „ 

y2 

A + a 1 

tn L_ 

i . a til 

1 A 

1 A yl 

tan Y* + 

1 - a 2 tan y* 


For a 2 - 0, tan y* ~ tan Y* + The greatest error would 

occur as y* approaches 0 or ir, where | Y' - Y*| approaches oy 
[ y * - y* | goes to zero as y* approaches ^ or 2 ' Similarly, 
for a 1 = 0, cot Y’ = cot y* - a 2 . The greatest error occurs as y* 
approaches ~ or where |y* - Y* I approaches o 2 . [Y r - Y*| 

approaches zero as y* approaches 0 or tt. We now have an indica- 
tion of the magnitude of error which is expected to be encountered. 

The stochastic aspects of this controller can be considered by 
the following simplified analysis: 
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Let y(t) = Asin (Wjt + <J>) + N(t) where N(t) is the component at 
frequency w of a white, Gaussian process with zero mean. 

y(t) = A ! cos Wjt + A^ 2 sin lO^t + N(t) 

Let Uj(t) = y(t) cos co^t 

= A , cos^unt + A - sin w.t cos ant + N(t) cos ant 
yl 1 y2 1 1 1 

= yA yl Cl + cos2co 1 t) + ~ ky 2 sin2u 1 t + N(t) cos a^t 
EfUjCt)] = E[jA yl (1 + cos2 Ul t) + jA y2 si^c^t] 

+ E [N (t) ]E [cosWjt] 

= jAyl (1 + cos2a) 1 t) + J A y2 Sin2(0 1 t 
Similarly, E[U 2 (t) = y (t) sina^t] = j A siirtt^t + \ \ 2 C 1 ' cosl^t) 


Thus, the expected values of the information signals are their deter- 
ministic components. The variance of the random component is: 


{ S v ,)2! - e[n 2 W] e[{^“;'} 2 ] 


1 n 2 

= 7 a N 


coi i 


where a., is the variance of N(t) 
N 


Thus the information signal y 
1 2 

— a N . This provides adequate 
aspects of the controller. 


has a variance of approximately 
information about the stochastic 
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5.3 Analog Simulation of Control System 

The circular perturbation system was simulated on the analog 
conputer. The system flow chart is shown in figure 66. An avail- 
able single output variable frequency sine wave generator provided 
the perturbation signal. The complement signal (cos w^t) was ob- 
tained by a simple approximate derivative circuit. It is recommended 
that the hard-wired unit be a single frequency oscillator providing 
both signals. Special purpose elements are given in figure 67. The 
plant was modeled as a quadratic cost function surface to enable com- 
parison of performance with other adaptive methods. Since phase 
response of the plant is equal in the two parameter dimensions, the 
two pure time delays were placed parallel to the plant. 

j = j (x 2 (t) + x 2 (t)). 

12 2 ^ 2 2 
= y (x^(t) + 2 Ax^ (t)cosoj^t + A cos^co^t + 2Ax^ (t) sinw^t + A~sin w^t) 

= j (x 2 + x 2 + A 2 + 2A (Xj(t)cosw 1 t + x 2 (t) sinojjt) ) 

Since z^ft) = ~ A^(t) forms the denominator for two division 
operations, we must guarantee that it never goes to zero to prevent 
circuit overloading. This can be accomplished by the introduction 
of a bias voltage. To minimize performance degradation, it is desir- 
able at activate the bias voltage only when z 3 ( t ) is small. Noting 
that such a bias voltage increases system damping, the following 
algorithm has been formulated: 
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z^(t) = z 3 (t) for z 3 (t) > B 

= m z 3 (t) + (l-m)B for z 3 (t) < B 
This is illustrated in figure 64. 



FIGURE 64 
Bias Voltage 

The slope m of the segment for z 3 C t ) < B > and B -> the switchiT1 S 
point, add two parameters which are to be manipulated for optimal 

system performance. 

There are two aspects to be considered in the formulation of a 
performance criterion. First is that of steady state error, e^. 

For the quadratic plant, e ss occurs when X^t) = X 2 (t) = 0. Then, 

, 4 » 2 - 

ss 2 

This non-zero steady state error is inherent in circular perturbation. 
The concept is illustrated in figure 65. 
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SPL 



Steady State Error Due to Circular Perturbation 

Let e = IP . Second is the convergence performance. Let IP- » 
ss 1 4 

/" ( e (t) _e ss ) 2 dt serve as a convergence performance where: 

e(t) = j t) + x^OO + A 2 + 2A( JCl (t)cosw 1 t + x 2 (t) sinWjt)) 

The inclusion of the sinusoidal components accounts for the performance 
degradation due to the perturbation of the plant. The total performance 
criteria can then be given by: 

IP = aj IP 2 + a 2 IP 2 

00 2 

= a l e ss + a 2 f o (eW - e ss } dt 

where Sj and a 2 are weighting factors. 

Optimization is to be accomplished by varying four parameters: 

A, the perturbation amplitude, K, the velocity gain, B, the switching 
level, and m, the slope of the damping bias voltage. Optimal values 
of these parameters were selected by forming a sample test array and 
observing general trends. 
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\ 


For K small enough so that the system comes to rest at the 

minimum, e was found to depend only on A, as predicted by 
ss 

equation 5.3-1: 



Calculated versus simulated steady state errors are shown in Table 12. 


TABLE 12 

Steady State Error 


A 

1 

e = -s- A 
ss 2 

e gs (simulated) 

i 

.01 

.00005 

1 - ™ ' 

.0003 

.05 

.00125 

.0012 

.10 

.0050 

.0049 


All quantities are given in machine units, the operational units of 
the analog computer. Note: the maximum is one machine unit which 

corresponds to ten volts. Thus, a machine unit corresponds to the 
normalized units on which this system is based. Steady state error 
shall be considered important and therefore weighted heavily. There- 
fore A will be minimized to the limit allowed by system noise. 
Values of K for which the system does not come to rest at the mini- 
mum will be considered unacceptable. 

To reduce the dimension of the test set, A was set to .01 
machine units. Tests were then conducted in real time over the 



CQ 
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sample space defined by K, B and M. The system was allowed to come 
to equilibrium at a selected initial input parameter setting, X(0) = 
(.8,. 8). Then the adaptive controller was set in operation by switch- 
ing the mode of the input parameter velocity integrators from initial 
conditions (IC) to operate (refer to Figure 66). The performance 
indicator network was activated simultaneously. A sample test set is 
given in Table 13. 

The parameters m = 0, 0 = .0050, and K = 10 were selected 
because of the observed performance, consistency, and stability 
characteristics. The index of performance weighting factors and 

a 2 in 

IP = a 1 IP 1 ♦ a 2 IP 2 

were not formally set. Tbe steady state error, IP^ was considered 
separately and given intuitive weighting far exceeding that of IP 2 - 
The results obtained for IP 2 are considered satisfactory despite 
this disproportinate weighting. 

Results for the selected parameters are shown in Figures 68 to 
70. A was maintained at .01 as a base for further testing. As 
shown in Figure 69, two-thirds reduction of the cost function J 
for the test initial conditions was obtained in .08 seconds. Steady 
state operation was obtained in .6 seconds. These are promising 
results for the application of automatic control to the SMS. 

Also subject to evaluation was the system response to a load 
change. This was simulated by changing the location of the minimum 
the (X r X 2 ) plane. So that comparative data could be obtained. 


on 



FIGURE 66 


Schematic for Analog Simulation of 
Circular Perturbation System 
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FIGURE 67 


Details of Special Elements 




Pure Time Helnv 









TABLE 13 


Sample Test Set 


Initial conditions: x^CO) = x 2 (0) - .8 machine units 


A = .01 

K 

IP 1 

IP 2 

3 

ii 

o 

1 

.0003 

.0967 

B = .0010 

5 

.0003 

.0254 

m = 0 

1 

.0003 

.0979 

B = .0050 

5 

.0003 

.0265 


10 

.0003 

.0138 

m = 0 

1 

.0003 

.1725 

B = .0100 

5 

.0003 

.0354 


10 

.0003 

.0153 

m = -10 

1 

.0003 

.0951 

B = .0010 

5 

.0003 

.0258 


10 

.0003 

.0139 

m * -10 

1 

.0005 

.1332 

B = .0050 

5 

.0005 

.0262 


10 

.0005 

.0139 

m = -10 

1 

.0005 

.1960 

B = .0100 

5 

.0005 

.0402 


10 

.0005 

.0154 
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FIGURE 68 

X(t) for _x(0) = [®] volts 


FIGURE 69 

X(t) for x(0) = fg] volts 




Scale = 4 volts/cm 



Scale = 4 volts/cm 




FIGURE 70 



Ordinate = J (SPL) 
Scale = 2 volts/cm 

Absicca = Time 


Scale *= .1 sec/cm 
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the load change was modeled as a step change from the minimum at 
X = (0,0) to a minimum location at X = (.8, .8). The parabolic 
surface contour was maintained. Since the final value for the load 
change test was the initial value for the "convergence" tests, a 
meaningful comparison of the two types of response could be obtained. 

Circuit overloading and instability resulted except for very 
small K and large B. Examination of the system indicates that the 
two dividers which yield sin y* and cos y* overload when the step 
change is made. This can be attributed to incompatible step change 
response of the parallel circuits providing the numerators and denomina 
tors of the two dividers. The circuit providing the denominators 
responds slower than the other two circuits (refer to Figure 66). 

This allows the numerator to exceed denominator, thus overloading 
the division circuits. 

Various corrective measures were tried, such as increasing the 
damping in the numerator circuits to obtain computable step response. 
But, acceptable performance to load change could not be obtained. 
Therefore, modification of the control scheme is necessary. 
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5.4 Modified Circular Perturbation Controller 

The following perturbation adaptive controller is proposed: 



FIGURE 71 

Modified Circular Perturbation Controller 
G (s) and G_(s) are pure time delays so that the phase shift across 

J. o 

{ G l*- s -*} and A. is exactly £2tt, £ = 0,1,2,... . G 4 represents servo 

GjCs) 

drives which will be utilized in the control of the physical plant. G 2 (s) 
is to be determined. 

The signal analysis is as follows: 
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J(t) * J^Ct) + A y Sin ( W 1 t + 


= + A ., cos ( 0,1 + A 0 sin w,t 


V 1 


For the outside loop, 


Y ? Ct) = (J dc (t) + A 1 cos (Ojt + A^ 2 sin o^t) cos o^t 


Y y (t) = J dc (t) cos Wjt + j A l (1 + cosZo^t) + j P ‘y2 sin2tL 'i t 


Y 7 (t) is operated on by G 2 (s). Assuming for now that G ? (s) = 
1, Y g (t) = Y ? (t) . Yg (t) is then integrated. 

v 8 (t) - K / Vt)dt 


For J, (t) slowly time varying, 
do 


yi 


Y 9 (t) = K t + 


J dc (t) \l 

sincOjt + sinZio^t 


to 


1 


'1 


A y ? 

~r~— COS2(0, t) . 
4(0. 1 


Note that integration attenuates the sinusoidal components by the 

A j 

amount of their frequency. K-^— t is the desired signal, so an 
improved signal would be obtained with G 2 (s) being a low pass 
filter. For our specific application, let 

G 2 (s) = G 4 (s) = •— j-^j- . 


A similar analysis holds for the inside loop. Again, let 
G 2 (s) = G 4 (s) = . 


For this system, the control gain in the system dimensions is 
proportional to the magnitude of the respective gradient components. 
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Thjs system more closely approximates the classical perturbation 
type gradient controller than the trigonometric form investigated 
in the previous section. The classical features will be investi- 
gated in section 5.5 to provide a basis for evaluation of system 
performance. 

The modified circular perturbation controller does not contain 
division circuits, so the circuit overloading due to load change 
experienced in the previous controller is not expected. Because of 
the different nature of this controller, it is difficult to predict 
its relative performance. 

This system was simulated on the analog computer. The flow 
chart is shown in Figure 72. For the quadratic plant the phase 
shift is zero so 

GjCs) - G 3 (s) = 1. 

A was set at .05 machine units to reduce noise effects. This was 
necessary because there is less signal filtering in this circuit 
than previously. K — 300 was found to optimize both convergence 
performance and load change performance. A sample data set is given 
in Table 14. The steady state error was .0009 due to the increased 
perturbation amplitude. Instability did not occur until K - 800, 
The index of performance shows that load change stability of 
the modified system was obtained at the expense of initial conver- 
gence response degradation. Two thirds reduction for the conver- 
gence testcase now requires .26 seconds as opposed to .08 seconds. 



l?3U 


FIGURE 72 

Schematic of Modified Perturbation System 
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TABLE 14 

Test Set for Modified Controller 

Cdnvergence Test 
A = .05 

Initial Parameter Setting :• ;x(0) = (0.8, 0.8) 
Optimal Parameter Setting : X* (T) = (0.0, 0.0) 


K 

IP 1 

IP 2 

100 

.0009 

.0679 

200 

.0009 

.0437 

300 

.0009 

.0366 

400 

.0009 

.0401 


Load Change Test 
A = .05 

Initial Optimal Parameter Setting : X_* (0) = (0,0, 0.0) 
Final Optimal Parameter Setting : ,X*(T) = (.8, .8) 


Functional 

Change : 

Step 

K 

IP 1 

^2 

100 

.0009 

.0847 

200 

.0009 

.0673 

300 

.0009 

.0670 

400 

.0009 

.0966 
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The same reduction for the load change test requires .40 seconds. 
But, this is still adequate system performance. Examples of system 
response are given in Figures 73 to 76. 



FIGURE 73 


Optimal Convergence Response 


FIGURE 74 


Optimal Convergence Response 
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X 


2 



Scale = 4 volts/cm 
x T (t=0) = [8,8)1 volts 
= [0,0] volts 




0 


Ordinate = SPL 

Scale = '2 volts/cm 
SPL (t=0) = .64 volts 

SPL (t-*») = .009 volts 

Absicca = Time 

Scale = .2 sec/cm 





FIGURE 76 


Optimal Load Change Response 


FIGURE 77 


Optimal Load Change Response 
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Scale = 4 volts/cm 
x T (t=0") = [0,0] volts 
x T (W) = [8,8] volts 



0 


Ordinate 
Scale 
SPL (t=0 + ) 
SPL (t-*») 
Absicca 
Scale 


SPL 

2 volts/cm 
.64 volts 
.009 volts 
Time 

.2 sec/ cm 
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5.5 Classical Perturbation Controller Features 

Eveleigh (9) provides the basic form of a classical single 
dimensional perturbation controller. Figure 77 shows the system 
with notation consistent with that of this chapter. 


X q( s ) 



FIGURE 77 

Single Dimensional Perturbation Adaptive Controller 


The plant is represented by" A(x) . A linear filter G^(s) represents 
the index of performance measurement process. G^(s) is generally 
a narrowband filter centered at frequency Experimental evidence 

(9) shows that a reasonable form for G^(s) is 


Gj(s) 



w 1 2 S 


2 

W 1 S + W 1 ' 


£[y 3 (t)] 


G 3 (s) 


is chosen to provide a phase shift at equal to that across 
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A(x), G m (s), and GjCs). The signal at Y 5 is approximately 
2 A 

± Ay sin^u^t = ± — if (1 - cos2Wjt). 

Since the desired signal is a > G,(s) is chosen as a low pass 

Y' 2 

filter. The integration provides a type-one adaptive loop, allowing 
the system to follow changes in the optimum parameter setting with 
no steady state error. 

Other than the concept of the perturbation circle, the narrow- 
band filter G^Cs) is the only basic difference between the classi- 
cal form and the controller developed in Section 5.3. This filter 
was incorporated into the controller so that performance data could 
be obtained. The filter circuitry was developed by the "M" method: 


Gj (s) 




+ to^S + C0 1 


2 " X(s) 


then. 


Y(s) 

ws 


tOjXCs) 


2 2 
S + (d^S + 


= M(S) 


2 2 

(s + to^s + toj ) M(s) = ==> 


2 

— M(s) = (j^XCs) - sM(s) - OJ.M(S) 
1 1 


and Y = to^sM(s) . 

The circuit can then be shown as in Figure 78. 
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FIGURE 78 

Implementation of Narrowband Filter G^(s) 

The high gains are due to u) = (60 Hz) (2 tt) = 376.5 rad/sec. 

The system of Section 5.4 was modified to include G^(s) and 
various low order low pass filters for G^(s) and G^(s) . This 
revised system is shown in Figure 79. The performance index IP^ 
was found to be an order of magnitude greater than that of the 
original system for all combinations. This occurred as a result 
of lightly damped system response due to G^(s), necessitating 
the use of low gains. Thus, this system modification is rejected. 



FIGURE 79 


Schematic for Analog Simulation of 
Modified Perturbation System with 
Bandpass Filter 
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Case IH 
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VI Discussion of Results 

The development and analysis of the wave equation for the 
Spinning Mode Synthesizer has verified the predicted optimal solu- 
tion. The optimal solution for the minimization of a ducted spinn- 
ing mode generated by a circular transducer array given by 


Z . 

03 


r . 

oj 


= axial position of j th transducer 
= Z q , j = 1,2, . . . ,n 
= radial distance of j** 1 transducer 


r o> j — 1 ,2 , . . . ,n 


9 0j 


. th 


= angular position of j transducer 

= “ 0-1), 3 = 1,2, ...,n 

Q . = source strength of j** 1 transducer 
' °3 


4 > • 


— , j = 1,2, ... ,n 

th 

= acoustic phase of j transducer 
= 0-1), 3 = 1,2, ...,n 


(1.3-1) 


(1.3-2) 


(1.3-3) 


(1.3-4) 


by a circular cancelling transducer array positioned at 


- 2^, j = 1,2, ... ,n 

r oj “ V j = 1 » 2 ’“*' n 

9 oj = (j_1) + b j “ 1,2> ' *‘ ,n 


(1.3-5) 

(1.3-6) 


is the predicted solution 


03 


= source strength of j ^ transducer 
= Q q , j = 1,2, . . . ,n 


(1.3-7) 
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a . = acoustic phase of j 
miT . „ . . , 

= — C2J-1) + Tt, J » 


th 

1,2 


transducer 
» • • • 


(1.3-S) 


Thie optimal solution has been shown to provide cancel lation for 
a finite range of m,p, and frequency which is governed by the set 
{M>. Table 11 is repeated for basis of discussion. 


TABLE 11 


{M(m o )} for the 4 Transducer Array Set for Cancellation 



{M} 

o 

II 

e 

{±4', ±12, ±20, ±28, . . .} 

0 


m = 1 
0 

{-3, 5, -11, 13, -19, 21, ...} 

m = 2 

{-2, 6, -10, 14, -18, 22, . . .} 

0 


3 

II 

w 

{- 1 , 7, -9, 15, -17, 23, ...} 

o 



Cancellation is effective until the lowest cutoff frequency of the 
lowest order circumferential mode in the set {M} is reached. Thus 
it was predicted that effective sound pressure level reduction could 
be obtained only over the frequency ranges shown in Table 15. 
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TABLE IS 

Frequency Range of Effective Cancellation 


f range 


V- = 

0 

II 

y = 2 

• • • 

m = 0 
0 

0 

to 

1378 

1378 to 1913 

- 

- 

m = 1 
o 

662 

to 

1511 

- 

- 

“ 

m = 2 
o 


- 


- 

- 

- 

m = 3 
0 


- 


- 

- 

- 









0 Hz = cutoff frequency of m = 0, y = 0 mode 

1378 Hz = cutoff frequency of m = 0, y = 1 mode 

1913 Hz = cutoff frequency of m = 4, y = 0 mode 

662 Hz = cutoff frequency of m = 1, 0 = 0 mode 

1511 Hz = cutoff frequency of m = 3, y = 0 mode 


Note: These results are valid for a 12 inch diameter duct. 
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A frequency margin of about 100 Hz to 200 Hz below the uppeT fre- 
quency limits is considered adequate to assure minimal degradation 
of sound pressure level reduction due to low decay rates in the 
axial direction of the limiting higher order mode. 

Satisfactory sound pressure level reduction has been obtained 

for the m = 1 circumferential mode below 1400 Hz, substantiating 
o 

the predicted frequency limit. Figure 21, which gives the results 
for a test case at 1560 Hz, provides evidence of the generation 
of an m * 3 mode above its cutoff frequency of 1511 Hz. The 
radial distribution of the spinning modes generated by the individ- 
ual arrays is characteristic of the waveform between the y - 0 
and y » 1 cutoff frequencies, having some characteristics of both 
distributions. But, the cancelled radial distribution is clearly 
that of a y = 0 spinning mode. This is strong evidence of the 
propagation of an m = ^3, y = 0 mode. The reduced level of this 
mode is expected due to energy loss to the cancelled m = 1, y = 0 
mode . 

Examination of the relative acoustic phases of the transducers 
gives further evidence to the generation of a higher order mode. 
Table 16 gives the acoustic phases for generation and optimal can- 
cellation of the m Q = 1 mode. 

TABLE 16 


Acoustic Phases for Optimal Cancellation of the m Q = 1 Mode 



, , ,, n 5 tt it 7 it 1 

Acoustic phase (rad) °~42’~4 7r 42 4 
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Evaluation of equation 1.3-4 for n = 8 and m = -3 yields the 
identical acoustic phases! Therefore, the combined arrays behave 
as an m = -3 generator. This phenomenon can be verified for all 
m which are elements of {M} given in Table 11. 

Sound pressure level reduction of the plane wave within the 
specified frequency range had been verified by test cases. The 
propagation of a higher order mode above 1800 - 1900 Hz is ill- 
ustrated by Figure 80. An m^ - 0, y = 1 mode was generated at 
2000 Hz, which is above the m = 4, y = 0 cutoff frequency at 
1913 Hz. The cancelled radial distribution is clearly that of a 
y ss o spinning mode. Thus, the m - ±4, y = 0 mode predicted by 
equation 3.3-F has been generated. 

A mode m = -2 has been predicted to propagate instead of a 
cancelled distribution when sound pressure level reduction of the 
m^ = 2 mode is attempted by the optimal solution. This was veri- 
fied and is illustrated in Figure 81. The cancelled distribution 
is an identical waveform 6 dB greater than the individual dis- 
tributions. This is a result of both arrays contributing to the 
excitation of the m = 2 mode. The acoustic phases of the trans- 
ducers are those given by equation 1.3-4 for m = -2 and n = 8, 
as expected. 

Excellent agreement has been obtained between the wave equation 
theory and experimental results for the range of m, y and frequency 
over which effective cancellation can be obtained. 





FIGURE 80 

SPL Versus Radial Distance for m = 0, p = 1 
Mode at 2000 Hz 
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FIGURE 81 


SPL Versus Radial Distance for 
m » 2 , y = 0 Mode at 1500 Hz 




RADIAL DISTANCE (in.) 
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The explanation of the limited range of acceptable cancellation 
performance does not challenge the test results of Chapter II. 

Variation of transducer response, if uncompensated, can produce 
substantial performance degradation. Because transducers are re- 
latively low impedance acoustic sources, circumferential variation 
of in-duct acoustic loading can also result in significant performance 
degradation. 

This performance degradation may prove to be insignificant for 
many applications. The reduction of the discrete tone peaks observed 
in a sound pressure level versus frequency plot for a fan-duct system 
to the flow noise level is perfectly satisfactory performance. Further 
reduction will have negligible effect since it will be masked by the 
flow noise. Thus cancelling array performance requirements depend on 
the specific application. 

If the performance degradation is at an acceptable level, the 
cancelling array can be governed strictly by the simple relations 
developed for a two dimensional cancellation system in Section 3.4. 

The reduction of the cancellation system dimensionality from 2n, 
where n is the number of transducers, to 2 greatly simplifies 
system optimization and thus is very desirable. Transducer quality 
specification and symmetry specification in the duct design should 
be considered in the overall design concept so that the desired 
cancellation can be achieved. 

Amplitude and phase compensation at each transducer can be 
employed to enable two dimensional control for cases of unsatisfactory 
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performance degradation. Fans are relatively high impendance acoustic 
sources and therefore would be effected little by small variations of 
the acoustic impedance of the duct in the circumferential direction 
at the plane of the source. Transducers, on the other hand, are 
relatively low acoustic impedance sources. Also, the variation of 
transducer performance can be expected to exceed the variation of 
source strength for the various generating mechanisms of a fan. Of 
significance are the individual rotor stator interactions and rotor 
blade pressure fields. So that the cancelling array generates a 
waveform of the proper symmetry, each transducer can be calibrated 
so that at its input would be correct for the desired acoustic output. 

The development of a measurement system to provide an index of 
performance remains a difficult problem. The sound pressure level 
inside the duct is to be minimized so that the effective perceived 
noise level (EPNdB) of the sound radiating from the duct is minimized. 
The form presented in Section 3.3, given below, is most easily applied 
in the presence of a single mode. 


IP = index of performance 


w 


- g. SPL (t) 

w i=l i i - 


where w = number of probes 

g. = weighting factor related to position of probe 
SPL. = sound pressure level at i position 


The probe positions and the corresponding weighting factors will 



depend upon the radial distributions of the mode to be reduced, 
assuming m constant, A pattern recognition technique can be 
applied to determine w and Consideration must be given 

to the "point cancellation" effect so that the index of perfor- 
mance produces the desired minimum. 

When more than one mode is present, pattern recognition should 
be utilized to isolate the modes which are to be reduced. For an 
actual fan, this can be difficult if there is significant energy in 
more than one mode for a given harmonic of the blade passage frequency. 

Adaptive control techniques have been found to be appropriate 
for maximizing cancellation of spinning modes of sound by the 
Spinning Mode Synthesizer. Experiment and analysis have verified 
the predicted optimal solution (see Section 1.3) over the range of 
m, p, and f for which effective cancellation can be achieved. 

Direct application of the optimal solution utilizing the input 
parameters to the source array would enable application of a simple, 
model based controller. But, such a design would not be consistent 
with the control system design objectives (see Section 3.1). Thus 
an unknown source strength and circumferential and axial phase angles 
have been assumed. The design process has been targeted toward final 
application by the use of adaptive control. 

The reduction of the system dimensionality to two, an amplitude 
parameter and a phase parameter, by assuming the general form of the 
optimal solution (equations 3.4-1 and 3.4-2) has greatly simplified 
the control problem. Control system design was based on the follow- 
ing features: 
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1) 2-dimensional plant 

2) unimodal cost function surface 

3) time varying minimum 

4) noisy cost function surface due to SPL measurement algorithm. 

A perturbation signal, gradient technique is well suited for the 
design criteria. The plant is a high frequency electrical - acoustical 
system, enabling high perturbation frequencies and thus adequate time 
response . 

Control system simplicity and reduced time response were obtained 
by the generation of a ’’perturbation circle". Plant output signal 
separation was obtained by multiplication by the two reference sinu- 
soids and low pass filtering (see controller algorithms in Sections 
5.2 and 5.3). The best overall system response was obtained with 
simple integral feedback (Section 5.4). Step change response on a 
normalized, centrally located, quadratic cost function surface satis- 
fied feasibility requirements. Two thirds reduction for the conver- 
gence test case with (X^O), X 2 (0)) = (.8, .8) required .26 seconds. 
A step load change moving the minimum from (0,0) to (.8, .8) 
resulted in a two thirds reduction in .40 seconds. On line response 
of the same magnitude can be expected if the transducer voltage and 
phase are electronically controlled. 
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VII Summary, Conclusions and Recommendations 
7 . 1 Summary 

The reduction of the discrete tones generated by the fan, com- 
pressor, and turbine sections of jet engines is essential for jet 
aircraft to meet present and proposed noise rules. The EPNdB rating, 
where pure tones are penalized, would respond proportionately to 
reductions in the discrete tone levels. Present reduction methods, 
such as the installation of acoustical treatment, have high cost, 
as well as severe weight, efficiency, and maintenance penalties. 

The reduction of the discrete tones, or spinning modes, by the 
cancellation effect due to the combination of two acoustic fields having 
equal amplitudes but opposite phase has been proposed by the Noise Con- 
trol Laboratory at The Pennsylvania State University. Cancellation 
would be maximized by the optimization of an amplitude parameter and 
a phase parameter of another spinning mode source. 

The Spinning Mode Synthesizer has provided the means for effective 
study of this noise reduction scheme. A modified SMS is the proposed 
source of the cancelling acoustic field in actual application. For the 
purposes of this feasibility study, the synthesizer consists of 
electrical-acoustical transducers located in an equally spaced circular 
array. Two arrays are obtained by assigning alternate transducers to 
each array. A specified spinning mode is generated by the "source" 
array. The identical mode with appropriate phase and amplitude for 
cancellation is then generated by the "cancelling" array. 
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Open loop control of the SMS has been found to be unsatisfactoiy 
for efficient parameter variation, rigorous evaluation of the synthe- 
sizer, and the introduction of complexities characteristic of jet 
engines. Adaptive control techniques have been found to be appropriate 
for automatic control to maximize cancellation. Design has been based 
on a 2-dimensional, unimodal, time varying, noisy cost function surface. 
A perturbation signal, gradient technique was found suited for the 
design criteria and developed. Analog simulation has shown favoiable 
convergence and load change response. 

7.2 Conclusions 

The feasibility of the cancellation concept has been demonstrated. 
The functional dependence of the range of circumferential modes, radial 
distributions, and frequency on the number of transducers and the 
specified circumferential mode has been established for single arrays 
and two interacting arrays. The generation of single, specified modes 
to specified frequencies can therefore be accomplished by determining 

the minimum number of transducers. 

A perturbation type, 2 dimensional controller for the closed loop 
control of the cancelling array of the Spinning Mode Synthesizer has 
been designed for eventual application to fan-duct systems. The good 
time response of the system obtained by analog simulation has shown 
the feasibility of application to time varying spinning mode sources. 
An example of such a source is the jet engine on takeoff and landing. 
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The development and analysis of the wave equation for the 
Spinning Mode Synthesizer has verified the predicted optimal solution 
for the minimization of the sound pressure level of a ducted spinning 
mode generated by a circular transducer array. The optimal solution 
has been shown to provide cancellation for a finite range of circum- 
ferential modes, radial distributions, and frequency. The frequency 
range of effective cancellation is limited by the set (M) of 
circumferential modes whose elements are given by 

m = m + (2k+l)n 
o 

where m^ = specified circumferential mode 
k = . . • > — 1 3 0 > 1, ... 

n = number of transducers in each array 

The propagation of the modes {M} due to excitation by the combined 
arrays masks the cancelled mode. Increasing n increases the range 
of effective cancellation. 

Substantial sound pressure level reduction was obtained for 
modes within the limits imposed by the sets of {M}'s. Examination 
of the radial distributions of cancelled modes above the frequency 
limits verified the generation of higher order modes. Sound pressure 
level reduction within the allowed range of frequencies was slightly 
hampered by effects due to the lack of symmetry of the test system. 
The combined effects of different response characteristics from 
transducer to transducer, a non-circular duct, and plate radiation 
resulted in degradation of the cancellation performance. 
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7.3 Recommendations 

It will be necessary to develop a measurement system to provide 
a cancellation index of performance to apply automatic control to a 
cancellation system. The weighted probe array discussed in Chapter VI 
is physically realizable and suitable for development. Pattern recog- 
nition techniques should be applied to determine probe positions and 
weighting factors. 

The controller which has been developed should be employed if 
extensive experimental research is considered. The development of 
the measurement system to provide a cancellation index of performance 
should improve open loop efficiency to an extent that it would be 
appropriate for small scale testing. The measurement system is 
considered essential foT productive testing at any scale. 

The application of cancellation to a ducted, single Stage fan by 
a modified Spinning Mode Synthesizer is the logical next step m the 
development of the cancellation system. The equations developed for 
a single transducer array can be used to specify the minimum number 
of transducers necessary to generate the appropriate modes to specified 
frequencies. Transducer quality and duct symmetry should be considered 
in the design of the system so that the desired cancellation can be 
achieved. It is recommended that a duct which is circular within 
1/100 of an inch be used for such research. These considerations will 
also enable two dimensional control of the transducer array without 
amplitude and phase compensation at the individual transducers. 
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APPENDIX 1 

General Solution of the Wave Equation for a Semi-Infinite, 
Constant Radius, Rigid-Walled, Cylindrical Duct 

The system is shown in Figure 82 


r 



FIGURE 82 

Cylindrical Coordinate System 

Solution is based on the following premises: 

(1) Forcing occurs only at the plane defining the closed end of 
the duct. 

(2) Forcing will be steady sinusoidal, represented by the complex 

exponential e such that 

_ r . - icot n . -iwt 

A cos tot = Re [A e ] ^ A e 
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Supplementary to premise 2, the following corollary is applied: 

(1) All other system variables can be represented by the form 
B (cj) cos (wt + = Re [B (i(o) e^ Wt ] -*■ B (iw) 

where 

— jV 

B (ia)) = B (a>) e 

The general complex solution is obtained in this Appendix. The 
real component is evaluated in Chapter III for a specific forcing 
function. 

The wave equation is: 


V 2 $ (r,0,z,t) = - - ? ■ ■ 


( 1 ) 


3t‘ 


where $ is the velocity potential and the V 2 operator for a cylindrical 
coordinate system is 


V 2 


3 _ 

3r 







assume (by the method of separation of variables) that a solution exists 
of the form 


$(r,e >Z ,t) = R00 © (9) Z(z) T(t) 

= ¥ (r,9,z) T (t) 

An exponential solution is desired for the time variable, so T(t) will 
be of the form 


T(t) 


A e^ Wt + ^ 



225 


Then substitution of $ into equation 1, yields 

1 (? TOO) 


V 2 ($ T (t) ) = 


at 2 


(V 2 $0 T(t) = - — $ T (t) 
c 2 


W 2 T 




The Helmholtz equation follows: 


V 2 $ + k 2 $ = 0 where k = — 


The expanded form is 

r L 3r 3r" r z 3e z 3z z 

With the assumed solution 

$ = R (r) 0 (6) Z (z) T(t) 
and the homogeneous boundary conditions 
1) 0 (9) = 0 (9 + 2tt) 


= 0 


2 ) 


30(9) _ 30(9 + 2tQ 


30 


39 


31 M 

1 St 


4) 

1 3r 


r=a 


r=b 


= 0 


= 0 


the solution components have been determined (for example, Morse (11)) 

J' (a k a b) 

R (r) = L [ J (k° r) - -2 jSL- Y (k° r) ] 

1 J 1 m my Y '(a k a b) m my 

m v mp 
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0 ( 9 ) 

Z (z) 

T (t) 


im0 

e + 


. a 
iy z 
» ' m U 


iwt , 
e + 


-im6 

e 


+ 



z 


-iojt 

e 


such that 

$ Jo Jo 0 (0) R Cr) Z (z) T (t) 


The boundedness condition is applied to Z (z) . For ^ 

is complex and positive. So that Z (°°) is finite, Cj = 0. The 
complete solution for the homogeneous boundary conditions is, therefore, 


L[L ff (k° r)] [A, e 1 ” 0 + A, e- 1M ] • 
m=0 y=0 L my v my 1 1 2 


-im0- 


• [C e 1Ym|jZ ] [D, e iwt ♦ D 2 


where 


J (ok° r) 

. a ,,a . r T r, a . m my v ,, a 

L my ^ k my r " [J m (k my r) " y « ( .a . m (k my 

m v my J 


and m = circumferential mode number, a non-negative integer 
y = radial mode number, a non-negative integer 
a = hub-tip ratio = 


w It. 

k° = eigenvalue of order y for the m mode having a hub-tip 
„ my 6 

ratio a 


v a = wave number in axial direction 
’my 



(k a 1 
v my' 


J , Y = Bessel functions of order m 
m’ m 
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For the case under consideration, the inner radius a is zero, 

j' (0) 

so o = 0 and — y 0. Therefore, 

to) 

m 


L (k r) = J (k r) 
mp v mp J m k my 

Propagation occurs only in the positive z direction, so the sign of 
iat is negative. The solution is then: 


$ = P J (k r) e 

m=0 p=0 mp m mp 


i(± + T m y Z - «t) 


This can be written as: 


$ - E ? [A e im9 + B e' im6 ] J (k r) e l( ^mp Z " wt ^ 
® " n^o p=0 L mp 6 mp e J m mp 


where A and B are amplitude coefficients to be determined, 
mp mp 

The particle velocity at the source provides an inhomogeneous 
boundary condition. The particle velocity, v, is given by 

3 13 3 

X = - v $ = - [£ r 97 + ie 7 36 + -z 3z] 

where a. are unit vectors. Also, the pressure is given by 

—l 

p = P !£ = i 03 P <D. 


3$ 

v = - a tt- 
— r — r dr 


^ 34 

^6 -^6 r 30 


Thus 



228 


dft 

v = - a -~— 
— z — z 3 z 


P = i 0) p 0 

The velocity at the source is given by: 

v g (r,6,z = 0) = f (r, 0) cos cot 

— z j z = 0 

which is represented in complex form by: 

- , Q s itot 
v = f (r,6) e 


Since v = v„ 
s — z 


v = - a rr 

—2 —Z o Z 


we have 

z = O’ 

3$ 


? f i V [A e im0 + B e~ in50 ] J (k r) e i( W 
m=0 y=0 'my L my my J m my 


- cot) 


and 1^ = 0 


? oo imfi -im0n — ■ f -■ s itot 

l i v FA e + B 1 J k „r e 

m=0 y=0 T my L my my 1 m r ' nn 


m my 


The amplitude of v is equal to 


f(r, e ) - - m l 0 yip 1 Y m y t A my e + B my e ^ J m ( ' k my r - 1 

- - (Jo ai " e p!o ^ * 

? v-ime- f Y B J fk r)} 

m=0 y=0 'my my m l my ;; 

For convenience, the amplitude coefficients and are defined to be 

C my ^ " y=0 Y my A my ^m ^ k my r ^ 

D my ^ ” y=0 Y my B my ^ k my r ^ ‘ 

Then, 


-im0- 


r 00 ^ _ ijr.0 r. f a -im0 \ 

f(r,6) - -i { Z 0 C fr) e * m I 0 D (r) o ) 


m=0 my 
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Applying orthogonality to the angular variable, 


/f .£(r,9,0) e ln6 de - -it % A C^We’ 1 ' 8 a^de 


* ; 2 ” E n D (r) a 1 " 6 e^de} 
o m=0 mp v 


r 2-TT „ „ , xm6 -in0 Aa 

f Z_ C (r) e e d0 
o m=0 mp 


for m i n 

2tt C (r) for m = n = 0 

mp v 


2tt C (r) for m = n f 0 
mp k 


0 


for n ^ m 


/ 2lT Z„ D (r) e‘ im0 e" in6 de = < 2tt D (r) for n = m = 0 
o m=0 mp^ J mp v 

0 for n = m ? 0 


Therefore, 


fCr,e,0) e—de . -i 2„ tC^tr) * e 0 D nu (r)} 


-im0 


= 1 for m = 0 


where e 


= 0 for m i 0 
Applying orthogonality again, 

/f f(r,0,O) e in0 d0 = -i iff j Q C my (r) e imB e +in6 d0 ♦ 


,2tt _ f v -im0 A in0, O ' l 

+/ Z„ D (r) e e d0> 
o m=0 mp v ' 


which results in: 


! 2 o nr.e.o) V r) * V r)} 


im0 
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For m = 0, 


for m 4 0, 


C (r) + D fr) = •=— f 2v f(r,0,O) d0. 
oy oy 2 tf o 


c (r) = I_/ 27r £(r,0,o) e" im9 de 

my v J 2tt o ^ ’ ’ ' 

D (r) = t- f 2v f(r,0,O) e ira0 d0. 
my^ J 2tt o 


Now, for m = 0, 


C fr") + D fr) = y [A + B ] J (k r) 
oy^ J oy 1 - ) y=0 'oy L oy oy J o oy 


h f ? f(r - e - 0) de 


Applying orthogonality to the radial variable yields 


/ b L rv [A + B ] J (k r) J (k r) dr 
o y=0 ’oy 1 oy oy J o oy o oq 


and 


= i_ / 2tt / b r f(r,8,0) J (k r) dr dG 
2 tt o o v ’ ’ o oq 


/ r J (k r) J Ck r) dr 
o o v oy 1 o ^ oq 


for n / q 


2k 


\ — [k 2 b 2 - 0] J 2 (k b) for y = q 
2 1 oy J o oy n 

"oy 


T b2j oXy b ) 


CRef. 12) 


Thus, 
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A ° W * B0li ‘ ’^Xp^op 


I 2 o 'o r f(r ' e - 0) V k oq r) 


For m 4 0, 


C my^ y=0 ^my A my J m^ c my r -* 

= Ijj: f 2 J f(r,0,O) e" ime d0 


Applying orthogonality again. 


^ o y=0 ^my A my J m ( ' k my r ' ) J m ( - k mq r ' 1 dl 

- H ! V ! l r f(r ' e ’ 0) e ' ime W 5 dr de 


r T J (k r) J (k r) dr 
o m my nr mq 


for y i q 


(Ref 


2k 


i_ [k^ b 2 - m 2 3 J k (k my l>) for P = q 

'my 


Thus, 


i k 2 
my 


rau TTY [k 2 b 2 -m 2 ] J 2 (kb) 

r my L my m my 


,2tt r b £f Q n . -im0 
/ / rf (r, 0,0) e 

o o 


J (k r) dr d0 
m ^ mu 


By the same procedure 


i k‘ 


B 


my 


tty [k 2 b 2 -m 2 ] J 2 (k b) 
my my m my 


r 2TT ,-b _ , _ n , im0 

/ / rf (r, 0,0) e 

o o 


• J (k r) dr d0 
nr my 


dr d6 


. 12 ) 


which completes the solution for a general harmonic source at z = 0. 
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The pressure distribution is given by: 
PCr,0,z,t) = i p w $ 


00 oo 


im6 . e " iTne i j r k jA 

my e J nr my 


■ 4 » » Jo v io tV“" * ] e m,J 


where A , B . and [A + B ] are given above. This can be rewritten as 
my* my’ oy oy 

P(r,e,z.t) - * J 0 p S 0 E mp [O w e'” 6 . % e’ 1 " 9 ] • 

• J (k r) 

nr my 


where 


= P (r,0,z,ito) e 


loot 


p 03 k 


my 


iry [k 2 b 2 -m 2 ] J 2 (k b) 
r my L my J m my 


G « f 2V f h rf (r, 0,0) J (k r) e~ m6 dr d0, m ^ 0 
my o o m my 


H my “ f V f l rf ^ 9 ’ 0) eim6 ^ d6, m ¥ ° 


G + H - / 2lT f b rf(r,0,O) J (k r) dr d0, m = 0 
oy oy o o o oy 


The above development was derived from analysis in unpublished 
notes by M. J. Os lac. Research Assistant at The Pennsylvania State 
University's Noise Control Laboratory. 

For convenience in early treatment in Chapter I, the symbols g + 
and g are employed where g + and g_ are as follows: 


• g = E G J (k r) e 
6 + my my m my 


i(y z - tot) 
'my 


g = E H J (k r) e 
- my my m my 


i (y z - tot) 
v 'my 



